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ABSTRACT 


This  thesis  reports  on  aethods  for  solving  traveling 
salesman  problems  with  time-window  constraints.  Two  types  of 
time  windows  are  considered:  hard  time  windows,  which  are 
inviolable,  and  soft  time  windows,  which  are  violable  at  a 
cost.  For  both  cases,  we  develop  several  heuristic  proce¬ 
dures,  including  some  that  are  based  on  Stewart's  [Ref.6] 
effective  heuristics  for  the  traveling  salesman  problem 
without  time-window  constraints.  In  addition,  we  develop 
exact  algorithms  for  each  case,  which  are  based  on  the 
state- space  relaxation  dynamic  programming  metnod  of 
Christo  tides,  Mingozzi,  and  Toth  £Bef.5],  Computational 
experience  is  reported  for  all  the  heuristics  and  algorithms 
we  develop. 
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I.  IHIBODOCTIOH 


A.  OVEBVIEW 

Consider  a  traveling  salesman  having  to  visit  n  cities 
or  customers.  He  starts  from  a  depot  and  needs  to  visit  each 
of  the  other  n-1  cities  only  once  and  then  return  to  the 
depot.  The  cost  of  traveling  between  any  pair  of  cities 
(expressed  in  terms  of  distance,  time  or  cost,  etc),  say 
from  city  i  tc  j,  is  given  as  c;j  in  a  cost  matrix  C.  The 
problem  is  to  design  a  tour  through  the  n  cities  that  would 
minimize  the  total  cost  of  the  tour.  This  is  known  as  the 
Traveling  Salesman  Problem  which  is  a  well-known  classical 
operations  research  problem. 

The  TSP  is  called  Euclidean  when  the  cities  that  must  be 
visited  all  lie  on  the  same  plane  and  the  cost  of  traveling 
between  any  pair  of  cities  is  the  Euclidean  distance  between 
them. 

The  TSP  is  an  NP-complete  problem  [ Bef .  1,  2].  All  known 
exact  solution  methods  have  a  rate  of  growth  of  the  computa¬ 
tion  time  which  is  exponential  in  n.  On  the  other  hand, 
heuristic  solution  methods  have  a  rate  of  growth  of  the 
computation  time  which  is  a  lew  order  polynomial  in  n  and 
have  been  experimentally  observed  to  perform  well.  For  this 
reason,  there  has  been  an  extensive  amount  of  research 
directed  at  TSP  heuristics. 

In  this  thesis  we  consider  adding  time  window 
constraints  to  the  TSP.  That  is,  if  t;  is  the  time  mat  the 
salesman  visits  city  i,  then  t^must  satisfy  1-L  <  t- <  u- , 
where  lp  and  u^  are  the  specified  lower  and  upper  bounds  of 
a  time  window.  This  problem  is  not  as  well  studied  as  the 
unconstrainted  TSP,  but  there  have  teen  a  few  approaches 
used  on  the  problem. 
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Psaraftis  [Ref.  3]  has  presented  a  dynamic  programming 
model  and  solution  procedure  for  two  dial-a-ride  problems, 
which  are  similar  to  time-window  constrained  TSPs.  Baker 
[Ref.  4]  has  presented  an  exact  algorithm  using  a  branch  and 
bcund  procedure  which  is  effective  for  very  small  n. 
Christofides  et  al.  [Ref.  5]  have  presented  a  dynamic 
programming  state  space  relaxation  procedure  to  compute 
bounding  information  within  a  branch  and  bound  algorithm. 

The  objective  of  this  study  is  to  develop  exact  and 
heuristic  algorithms  which  will  provide  an  optimal  or  near 
optimal  tour  that  visits  each  city  in  its  given  time 
window.  He  are  given  a  depot  location,  a  set  of  x,y 
co-ordinates  for  n  cities  and  a  set  of  time  windows  . 

A  common  application  of  the  TSP  is  in  vehicle  routing 
problems.  a  set  of  customer  orders  must  be  partitioned 
among  several  vehicles.  Given  a  partition,  the  problem  tnen 
decomposes  into  one  TSP  for  each  vehicle.  Because  of  this 
prospective  application  and  in  deference  to  the  difficulty 
of  large  TSPs  (with  or  without  time  constraints),  we  confine 
our  research  and  computation  to  small-scaled  problems  (less 
than  3  0  nodes). 

We  consider  two  different  kinds  of  time  window 
constraints;  hard  time  windows  and  soft  time  windows.  Hard 
windows  cannot  be  violated.  Soft  windows  can  oe  violated, 
but  a  penalty  cost  must  be  paid  when  they  are.  The  penalties 
can  be  defined  individually  for  each  customer,  and  they  can 
differ  for  early  and  late  arrivals.  Generally,  the  penalty 
for  arriving  before  the  lower  time  window  bound  is  much  less 
than  tne  penalty  for  arriving  after  the  upper  bound.  In 
Chapter  III,  we  present  tne  hard  time  window  approach  and  in 
Chapter  IV,  we  present  the  soft  time  window  approach. 

We  developed  several  Fortran  programs  for  solving  the 
TSP  and  time-constrained  TSP.  For  the  TSP,  we  use  Stewart’s 
[Ref.  61  recent  heuristics,  CCCO  and  CCAO.  For  the  time 
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constrained  ISP  problems,  we  develop  some  new  heuristics, 
some  of  which  are  modification  of  Stewart's  heuristics  for 


the  unconstrained  problem.  He  also  developed  exact  algo¬ 
rithms  for  both  hard  and  soft  windows  using  Christofides  et 
al. 's  £Ref.  5]  method  of  state-space-relaxation  dynamic 
programming  and  branch  and  bound.  This  is  described  in 
Chapters  III  and  IV. 

Finally,  we  describe  a  hybrid  of  the  heuristic  and  exact 
programs.  The  hybrid  uses  the  overall  structure  of  the 
exact  program,  but  the  upper  bounds  are  obtained  with  the 
heuristic.  Ibis  is  discussed  in  Chapters  III  and  IV. 

B.  THE  TRAVEIIHG  SALESMAH  PROBLEM 

A  tour  is  a  chain  which  passes  through  all  the  n  nodes 
and  in  which  the  first  and  the  last  nodes  coincide.  A  tour 
is  also  known  as  a  Hamiltonian  cycle. 

Let  a  tour  be  denoted  by  t  =  (i  ,  i  ,...,  i  ,  i  )  and 

12  n  1 

the  cost  of  this  tour  be 
n- 1 

C  (t)  =  c  ♦  c  . 

■*=- — >  XX  X  1 

j  =  1  j  j  +  1  n  1 

Here  <  i  ,  i  ,...,i  )  is  a  permutation  of  the  integers 
1  2  n 

from  1  to  n,  giving  the  order  in  which  the  cities  are 
visited . 

The  Traveling  Salesman  Problem  can  be  defined  as 
follows.  Given  a  graph  G  =  {  N,A  }  composed  of  a  set  of 
nodes  N  ,  a  set  of  arcs  A  connecting  these  nodes,  and  a  cost 
(distance  )  c$j  associated  with  each  arc  (i,  j)  in  A.  The 
TSP  is  the  problem  of  finding  the  minimum  cost  tour  of  the 
nodes  in  N.  The  following  mathematical  formulation  of  the 
TSP  is  from  Stewart  £Ref.  6  "| . 
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where 


if  arc  (i,  j)  is  on  the  tour 
otherwise 


y  are  continuous  variables  that  force  the  final 
ij 

solution  to  be  on  the  tour 

(  i.e.  include  every  node  in  the  same  route  ) 
u  is  a  number  >  n- 1  ,  and 

n  is  the  number  of  nodes  in  the  set  N. 


The  constraints  (1.2)  and  (1.3)  ensure  that  each  node 
will  be  visited  exactly  once,  while  constraints  (1.4)  and 
(1.5)  force  the  final  solution  to  be  a  single  tour  that 
starts  and  ends  at  node  1  (depot)  .  This  formulation  is  not 
directly  used  in  our  TSP  programs,  but  is  of  interest  in  a 
general  discussion  of  the  problem. 


C.  TSP  WITH  II HE  WIHDCW  CONSTRAINTS 


The  time-constrained  Traveling  Salesman  Problem  is  a 
variation  of  the  TSP  that  includes  time  window  constraints 
on  the  time  to  visit  some  of  the  cities.  The  hard  time- 

constrained  TSP  is  to  find  the  minimum  cost  tour  subject  to 

visiting  each  city  within  its  time  window. 

For  the  time-constrained  TSP  model,  we  define  a  contin¬ 
uous  nonnegative  variable,  t-L  ,  to  be  the  time  that  the 
salesman  visits  city  i.  Since  the  salesman  must  return  to 
city  1  (depot)  at  the  end  of  the  tour,  the  formulation 

includes  an  additional  variable,  tn+J  ,  the  total  time 

required  to  complete  the  tour. 

He  assume  that  a  complete,  symmetric,  nonnegative 
distance  matrix,  |Cij  1,  is  known  and  that  time  is  a  scalar 
transformation  cf  distance.  Thus  time  and  distance  may  be 
used  interchangeably. 

The  following  mathematical  formulation  of  the  TSP  with 
time  constraints  is  from  Baker  £Hef.  4]. 


MIN  t 

n  +  1 

t 

1 

(1-6) 

S.T  t 

i 

t  > 
1 

c 

1i 

i  =  2,3,. 

..  ,n 

(1.9) 

1  t  - 
r 

t  1  > 
3 

c 

ij 

j  =  3,4,.. 

2  <  i  <  j 

.  ,n 

(1.10) 

t  - 
n  +  1 

t  > 

i 

c 

1i 

i  =  2,3,.. 

.  ,n 

(1.11) 

t  > 

i 

0 

i  —  1,2,... 

,n+ 1 

(1.12) 

1  < 
i 

t  <  u 
i 

i 

i  =  2,3,.. 

.  ,n 

(1-13) 

where  t  =  the  time  that  the  salesman  visits  city  i 
i 

U  |  =  the  absolute  value  of  x 

c  =  the  shortest  time  required  to  travel  from 

ij 

city  i  to  city  j 
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1  =  the  lower  bound  on  the  time  window  for  the 

i 

salesman  to  visit  city  i 

by  assumption  all  1  >  0 

i 

u  =  the  upper  bound  on  the  time  window  for  the 
i 

salesman  to  visit  city  l 

u  >  1  ,  for  all  i 

i  i 

The  constraints  (1.9)  through  (1.12)  ensure  a  nonnega¬ 
tive  arrival  time  at  city  i,  ,  be  obtained  for  each  city 
(node  2  through  node  n)  on  the  tour.  The  constraint  (1.9) 
guarantees  that  t;  ,  the  time  that  the  salesman  leaves  the 
node  1  (depot)  will  be  the  smallest  t  value.  The  absolute 
value  constraint  (1.10)  ensures  that  the  arrival  times  of 
any  two  city  differ  by  amount  of  time  sufficient  to  allow 
the  salesman  to  travel  between  the  two  cities.  The 
constraint  (1.11)  guarantees  that  t,*,  ,  the  time  the  salesman 
returns  to  the  depot,  will  be  the  largest  t  value.  The 
inequalities  (1.12)  and  (1.13)  are  nonnegativity  and  the 
time  window  constraints  respectively. 

Unfortunately,  Baker’s  proposed  model  for  the  time- 
constrained  ISP  is  very  difficult  to  solve,  because 
constraint  (1.10)  is  nonconvex.  Therefore,  we  will  not  use 
this  formulation  directly  in  cur  program- 

Figure  1.1  illustrates  the  nonconvexity  of  constraints 

(1.10)  for  one  i,j  pair,  the  example  |  t3  -  t  |  >  5.  The 

feasible  region  for  this  constraint  is  the  union  of  two 

disjoint  sets.  Taken  all  together,  constraints  (1.10) 
m 

define  2  disjoint  sets  sets  where  m  =  (n- 1)  (n- 2) /2, 

which  are  very  difficult  to  work  with. 

We  can  see  that  the  time-constrained  TSP  is  very 
different  from  ISP,  even  in  formulation. 


II.  HEURISTIC  TSP  SOLUTION 


A.  OVEBYIEH 


Many  heuristic  procedures  have  been  developed  for 
solving  TSP.  Our  purposes  in  this  Chapter  are  to  examine 
some  of  the  well-known  heuristics,  to  review  Stewart's  [Ref. 
6]  recent  heuristic,  and  to  compare  these  approximate  tech¬ 
niques  cn  the  basis  of  efficiency  and  accuracy  on  a  small 
number  of  examples. 

In  general,  heuristic  procedures  are  categorized  by 
three  broad  classes:  tour  construction  procedures,  tour 
improvement  procedures,  and  composite  procedures  [Ref.  7]. 
Tour  construction  procedures  start  with  a  single  node  and 
successively  add  nodes  till  a  tour  is  built.  Tour  improve¬ 
ment  procedures  attempt  to  find  a  Jtter  tour  given  an 
initial  tour.  Composite  procedures  construct  a  starting  tour 
from  one  of  the  tour  construction  procedures  and  then  try  to 
find  a  tetter  tour  using  one  cr  more  of  the  tour  improvement 
procedures. 

1-  lour  Construction  Procedures 

There  are  many  methods  available  for  constructing  an 
initial  tour.  Procedures  which  have  been  generally  used  are 
given  below. 

a.  Nearest  Neighbor  (  Rosenkrantz  et  al.  [Ref.  8]  ) 

Step  1.  Start  with  any  node  as  the  beginning  of 
a  subtour. 

Step  2.  Find  the  node  closest  to  the  last  node 
added  to  the  subtour.  Add  this  node  to 
the  current  subtour. 


Step  3.  Repeat  step  2  until  all  nodes  are  contained 
in  the  tour.  Then,  join  the  first  and  last 
node. 

b.  Clarke  and  Wright  Savings  {  Clarke  and  Wright 

[Ref.  9]  ) 

Step  1.  Select  any  node  as  the  central  depot 
which  we  denote  as  node  1. 

Step  2.  Compute  savings  s  =  c  ♦  c  c 

ij  1i  1j  ij 
for  i,j  =  2,3,...,n.  i#j 

Step  3.  Order  the  savings  from  largest  to  smallest. 

Step  4.  Starting  with  the  largest  savings  on  the 
list,  subtours  are  assembled  such  that  the 
next  node  added  has  the  largest  remaining 
savings  -  provided  that  a  constraint  is  not 
violated.  Once  a  pair  of  nodes  i  and  j  have 
been  linked,  they  remained  linked  . 

Repeat  until  all  nodes  have  been  assigned. 

Here,  the  quantity  s.  _  is  the  amount  of  travel 
saved  if  node  j  is  visited  directly  after  i,  as  opposed  to 
having  separate  trips  from  the  depot  to  nodes  x  and  j. 
Figure  2.1  demonstrates  the  procedure  for  two  nodes  i  and  j. 

c.  Insertion  procedures  (  Rosenkrantz  et  al.  [Ref. 

8]  ) 

An  insertion  algorithm  constructs  a  feasible 
tour  by  successively  adding  one  node  to  an  existing  subtour. 
This  procedure  takes  a  subtour  of  k  nodes  at  iteration  k  and 
attempts  to  determine  which  node  not  in  the  subtour  should 
join  the  subtour  next  {the  Selection  step)  .  Then  it  deter¬ 
mines  where  in  the  subtour  it  should  be  inserted  (Insertion 


Figure  2,1  Concept  of  the  Clarke  -  Bright 
Savings  Beuristic. 

step)-  Stewart  [  Hef .  6]  presented  the  following  general 

algorithmic  structure- 

step  1.  (Initial  Sufctour) 

Obtain  a  IS P  tour  for  a  subset  of  the  nodes 
N'c  N  in  G. 

Step  2.  (Selection  Step) 

Find  a  node  k€  N-N'  to  be  added  to  the 
existing  subtour. 

Step  3.  (Insertion  Step) 

Choose  an  arc(i,j)  in  the  subtour  on  N*. 
Insert  node  k  between  i  and  j  and  add 
k  to 

Step  4.  If  N  =  N ' ,  then  stop- 

(We  have  a  Hamiltonian  cycle). 

Otherwise,  return  to  step  2. 
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There  are  many  variations  on  this  algorithmic 
structure  depending  on  the  procedures  chosen  for  executing 
steps  1 ,2  and  3 . 

Wiorkowski  and  McElvain  [Ref.  10],  Or  [fief.  11], 
Stewart  [Ref.  12]  and  Norback  and  Love  [fief.  13]  all  present 
insertion  algorithms  that  use  the  convex  hull  of  the  set  of 
nodes  N  for  the  initial  subset  N • .  Nemhauser  and  Hardgrave 
[fief-  14]  have  shown  that  there  exists  an  optimal  tour  for 
every  Euclidean  TSfi  in  which  the  relative  order  of  the  nodes 
on  the  boundary  of  the  convex  hull  is  preserved.  This  means 
that  the  optimal  tour  visits  nodes  on  the  boundary  of  the 
convex  hull  in  the  same  order  as  if  the  boundary  itself  were 
followed. 

Further  justification  for  the  use  of  the  convex 
hull  for  the  initial  subtour  is  shown  empirically  by 
Stewart's  [Ref.  6]  computational  experiment.  He  compared 
several  insertion  heuristics  both  with  and  without  the 
convex  hull  as  the  starting  solution.  The  results  show  that 
all  the  insertion  algorithms  are  improved  by  the  use  of  the 
convex  hull.  Some  are  improved  substantially,  others  only 
moderat  ely. 

Many  criteria  have  been  suggested  for  the  selec¬ 
tion  of  the  node  to  be  inserted  in  an  insertion  procedure. 

( 1)  nearest  Neighbor  (  Rosenkrantz  et  al-  [  Ref. 

8]  ).  Choose  the  node  k  that  is  nearest  a  node  in  the 
current  tour.  I.  e.  ,  find  k  =  argmin  c  s.t.  j£N-N', 
i«i*.  j  ij 

(2)  Cheapest  Insertion  (  Rosenkrantz  et  al. 
F  Ref .  8]  ).  Choose  the  node  k  that  may  be  inserted  at 
minimal  increased  cost.  X.e.,  find 

k  =  argnin  c  +  c  -  c  s.  t.  m  €  N-N* ,  i, j  e  N' . 


(3)  Farthest  Insertion  (  fiosenkrantz  et  al. 
r  Ref.  8]).  Choose  the  node  k  that  is  farthest  from  a  node 
current  subtoor.  I.e. ,  find  k  =  argmax  c  s.t.  j  €  N-N  • , 


f  Ref.  8]  ) 


(4)  Arbitrary  Insertion  (  Rosenkrant z  et  al 
Choose  node  k  randomly  from  among  N-N  * . 

(5)  Ratio  I  nsertion  (  Stewart  [Ref.  J2  ]  ). 
Choose  the  node  k  such  that  the  proportional  increase  in 
cost  is  minimal.  I.e.,  find  k  = 
s.t-  m  €  N-N',  i, j € N* . 


argmin  (c  ♦  c 
m  im  m  j 


)  /  c. 


(6)  Perpendicular  Distance  (  Wiorkowski  and 
ttcElvam  [Ref.  JO]  ).  Choose  the  node  k  that  is  closest  to 
an  arc  in  the  current  subtour. 

(7)  Ratio  Times  Distance  (  Or  [Ref.  J  J  ]  ). 

Choose  the  node  k  such  that  the  product  of  ratio  and 

distance  is  minimized.  I.  e.  ,  find 

k  =  argmin  ((c  +  c  )/c  )x(c  ♦  c  -c  ) 

m  im  mj  ij  im  mj  ij 


S.t.  m  €  N-N*  ,  i ,  j  €  N  * . 


(8)  Greatest  Angle  (  Norback  and  Love  [Ref. 
131  ).  Choose  the  node  k  and  arc  i,  j  such  that  the  angle 
formed  by  the  two  arcs  (i,k)  and  (k,j)  is  a  maximum.  I.e., 
find  k  =  argmax  angle{  arc(i,m),  arc(m,j)  }  s.t.  m  €  N-N* , 

ID 

i, j  €  N •  . 


The  insertion  criteria  that  have  been  used 
fall  into  two  categories.  [Ref.  6] 

1.  Cheapest  Insertion 


Insert  the  node  kfe  N-N'  between  those  two 
that  minimize  the  guantity 


connected  nodes  i,j  €  N ' 


2.  Identical  Insertion  and  Selection 


step. 


Do  selection  and  insertion  in  the  same 


2.  lour  improvement  Procedures 

The  best  known  procedures  of  this  type  for  the  TSP 
are  the  branch  exchange  heuristics  [Ref.  7].  These  branch 
exchange  heuristics  work  as  follows. 

Step  1.  Find  an  initial  tour.  This  tour  may  be  chosen 
randomly  from  the  set  of  all  possible  tours, 
cr  it  may  be  generated  by  one  of  the  tour 
building  procedures  above. 

Step  2.  Improve  the  tour  using  one  of  the  branch 
exchange  heuristics. 

Step  3.  Continue  step  2,  until  no  additional 
improvement  can  be  made. 


For  a  given  k,  we  define  a  k-change  or  a  tour  as 
consisting  of  the  deletion  of  k  branches  in  a  tour  and  their 
replacement  by  k  other  branches  to  form  a  new  tour.  A  tour 
is  k-opt  if  it  is  not  possible  to  improve  the  tour  via  a 
k-change.  In  general  the  larger  the  value  of  k,  the  more 
likely  it  is  that  a  k-opt  solution  will  be  optimal. 
Unfortunately,  the  number  of  operations  necessary  to  test 
‘all  k  exchange  is  proportional  to  nK  ,  where  n  is  the  number 
of  nodes  in  the  TSP.  Due  to  this  complexity,  values  of  k  =  2 
and  k  =  3  ace  most  commonly  used  [  Ref.  7].  The  2-opt  and 
3-opt  heuristics  were  introduced  by  Lin  [Ref.  15]  and  the 
k-opt  procedure,  for  k>3  was  presented  by  Lin  and  Kernighan 


[Ref.  16]  . 

Or  [Ref.  11  ]  has  designed  a  modified 
considers  only  a  small  percentage  of  3-branch 
This  modified  3-opt  called  oropt  by  Stewart 


3-opt  that 
exchanges. 
[Ref.  6] 
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considers  only  those  branch  exchanges  which  are  composed  of 
a  string  of  one,  two,  or  three  adjacent  nodes  being  inserted 
between  two  ether  nodes  in  the  current  tour.  By  limiting  the 
number  of  exchanges  that  are  considered  in  this  way,  Oropt 
requires  many  fewer  calculations  than  a  full  3-opt. 

Stewart  [fief.  6  ]  made  an  experiment  of  the  convex 
hull,  cheapest  angle  insertion  algorithm  (CCA)  which  will  be 
discussed  in  the  next  section  as  a  stand-alone  algorithm  and 
with  each  of  the  three  post-  processors .  The  algorithms  are 
designated  CCA,  CCA2 ,  CCAO,  and  CC A3  for  the  convex  hull 
cheapest  insertion  stand-alone,  with  2-opt,  with  Oropt  and 
with  3-cpt  respectively.  He  drew  two  conclusions  from  his 
computational  results.  First,  the  3-opt  requires  substan¬ 
tially  more  time  than  either  the  2-opt  or  tae  Oropt. 
Second,  the  2-opt  is  dominated  by  the  Oropt  and  the  3-opt  in 
quality  of  solution. 

In  computation  time,  Oropt  only  looks  at 
2 

approximately  3n  of  the  n  possible  3-opt  exchanges  on 
eaca  pass.  There  are  n  ways  to  select  the  first  branch, 
times  3  ways  to  select  the  second  branch,  and  n-2  ways  to 
select  the  third  branch. 

This  accounts  for  the  fairly  close  times  for  the 
2-opt  and  Oropt.  The  quality  of  CCAO  solutions  dominate 
CCA2  solutions.  On  the  other  hand,  there  is  little  cr  no 
difference  between  the  Oropt  and  3-opt  in  terms  of  solution 
quality  . 

Stewart's  main  conclusion  fror  the  above  experiment 
is  that  the  Orcpt  performs  as  well  as  a  3-opt  in  a  small 
percentage  of  the  computer  time  required  by  a  3-opt,  and  it 
should  he  preferred  tc  both  the  2-opt  and  the  3-opt  for 
Euclidean  TSP's. 
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3.  Composite  Proced are 

The  basic  composite  procedure  is  a  combination  of 
the  tour  construction  and  branch  exchange  procedures.  It  is 
obtained  by  appending  a  branch  exchange  procedure  to  the 
tour  construction  algorithm  as  a  post- processor.  The  proce¬ 
dure  can  be  stated  as  follows  [Ref.  17]. 

Step  1.  obtain  an  initial  tour  using  one  of  the 
tour  construction  procedures. 

Step  2.  Apply  a  branch  exchange  procedure  to  the 
solution  produced  by  the  step  1. 

Stop  wnen  no  further  improvement  can  be  made. 

The  composite  procedure  is  relatively  fast  computa¬ 
tionally  and  gives  good  results  [Ref.  18]. 

E.  CC&C 

1 .  Algorithm 

The  CCAO  algorithm  designed  by  Stewart  [fief.  6]  uses 
the  convex  hull  of  the  nodes  in  N  for  its  initial  subtour. 
Then  it  inserts  the  nodes  not  currently  in  the  subtour  where 
they  may  be  inserted  most  cheaply  (the  Cheapest  Insertion 
criterion).  it  selects  the  node  k  to  be  inserted  at  each 
iteration  according  to  how  large  an  angle  is  formed  by  the 
two  arcs  that  must  be  added  to  the  current  subtour 
(Selection  criterion)  in  order  to  insert  k.  Finally  it  uses 
an  oropt  to  make  local  improvement  on  the  tour  constructed 
in  the  first  stage.  CCAO  means  Convex  Hull,  Cheapest 
Insertion,  Angle  Selection,  Oropt. 

Algorithm  :  CCAO 

Input  :  Number  cf  nodes,  x  and  y  co-ordinates  of  all 
nodes. 

Output:  Ordered  list  of  tour,  total  cost. 
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Step  1  .  (Initial  Sub  tour) 

Find  the  convex  hull  of  the  set  of  nodes  N. 
Call  the  set  of  nodes  on  the  boundary  N*. 
Let  the  initial  subtour  be  the  nodes  of  N' 
in  the  same  order  as  they  appear  on  the 
convex  hull. 

Step  2  .  (Cheapest  Insertion) 

For  each  node  m  6  N-N* ,  find 

(i  ,  j  )  =  argmin  c  ♦  c  -  c 
m  m  i,j  im  mj  ij 

s.t.  i,  j  €  N1,  i,  j  :  connected. 


Step  3  .  (Greatest  Angle  Selection) 

For  the  next  insertion,  select  the  node 

that  maximizes  the  angle  between  the  arcs 

(i  ,  m)  and  (m,j  )  over  all  m  6  N-N'. 
m  m 

I-e,  find  k  =  argmax  angle  (  (i  ,m),(a,j  )  } 

m  mm 


and  j 
k  k 


and  add 


s.t.  m  €  N-N* . 

Insert  k  between  i 
k  to  S', 
step  4  :  If  N'  =  N,  go  to  step  5. 

Otherwise  return  to  step  2. 
step  5  :  Apply  an  Oropt  to  the  current  tour.  Stop 
when  no  further  improvements  can  be  found. 

End  of  algorithm  CCAO 

2.  Exa  m  nle 

Figures  2.2  -  2.4  illustrate  the  above  algorithm  on 
the  TSP  defined  as  test  problem  [  1  ]  in  Appendix  A.  First 
the  convex  hull  is  generated  for  an  initial  starting 
subtour.  This  subtour  consists  of  nodes  2,13,12,14,5,15,7,4. 
A  solid  line  marks  the  boundary  of  the  convex  hull  in  Figure 


Pigure  2.2  Initial  Subtour  and  Insertion. 

In  step  2,  each  of  the  interior  nodes  (1,3,6,6,9,10,11,16) 
is  associated  with  a  pair  of  connected  nodes  on  the  initial 
subtour  (the  dashed  lines  in  Figure  2.2).  In  step  3,  the 
dashed  lanes  that  form  the  greatest  angle  (closest  to  180°) 
identify  the  node  to  be  inserted  (node  10  in  this  example). 

Figure  2.3  shows  the  problem-  after  the  first  three 
insertions (  ncde  10,  node  1  and  node  S  in  that  order). 
Notice  that  some  nodes  not  in  the  subtour  are  associated 


Figure  2.3  Intermediate  Subtoar  and  Insertions. 

with  new  node  pairs.  Figure  2.4  shows  the  final  tour  for 
stage  one.  This  tour  is  now  passed  to  an  Oropt  post¬ 
processor.  In  this  case  the  tour  from  stage  one  appears 
from  inspection  to  be  optimal,  and  oropt  will  find  no 
improve  scent . 


Figure  2.4  Final  Tour  of  CCAO. 

3.  Computational  Results 

In  addition  to  CCAO,  CCCO  (Convex,  Cheapest, 
Cheapest,  Oropt)  has  been  coded  for  the  purpose  of  compar¬ 
ison.  The  only  difference  between  CCAO  and  CCCO  is  that 
CCCO  uses  the  cheapest  selection  criterion  instead  of 
greatest  angle  of  CCAO- 

Xe  used  Sedgewick’s  £ Eef .  19]  package  wrapping  algo¬ 
rithm  for  finding  the  convex  hull  (initial  subtour). 


Starting  with  some  point  (called  the  anchor)  that  is  guaran¬ 
teed  to  be  on  the  convex  hull  (say  the  one  with  the  smallest 
y  co-ordinate),  take  a  horizontal  ray  in  the  positive  direc¬ 
tion  and  sweep  it  upward  until  hitting  another  point-  This 
point  is  on  the  hull.  Then  start  at  that  point  and  continue 
sweeping  until  hitting  another  point,  etc.  The  package  is 
completely  wrapped  when  the  first  point  is  included  again. 
The  following  algorithm  finds  the  convex  hull  of  an  array 
1(1,  ...,n)  of  nodes,  the  node  L(n+1)  is  used  as  a  sentinel, 
that  is,  a  copy  of  the  first  node  which  is  used  to  signal 
completion  of  the  procedure.  The  variable  NH  is  maintained 
as  the  number  of  nodes  so  far  included  on  the  hull. 

Algorithm  :  Package  Irapping 

Input  :  Number  of  nodes,  x  and  y  co-ordinates  of  all 
nodes. 

Cutputr  Ordered  list  of  convex  hull  and  number  of 
nodes  included  on  the  convex  hull. 

Step  1  .  (Initialization) 

find  and  duplicate  anchor.  I.e.,  find 

NMIN  =  argmin  y  s.  t.  i  €  N  and  set 

i 

NH  =  0,  L(n  +  1)  =  L  (NMIN)  . 

Step  2  :  (Swap  nodes  NH  and  NMIN). 

Put  last  node  found  into  the  hull  by 
exchanging  it  with  the  NHth  node. 

NH  =  NH  ♦  1. 

TEfiP  =  L  (NH)  . 

L (NH)  =  L  (NM  IN)  . 

L(NMIN)  =  TEME. 

Step  3  :  (Compute  angle) 

Compute  the  angle  from  the  horizontal  made 


by  the  line  between  L(HH)  and  each  of  the 
nodes  not  yet  included  on  the  hull. 

Step  4  :  (Find  next  hull  node) 

Find  the  node  whose  angle  is  smallest  among 
those  with  angles  bigger  than  the  current 
value  of  the  'sweep*  angle  (the  angle  from 
horizontal  to  the  line  between  L(NH-1) 
and  L ( NH) ) . 

Step  4  :  Stop  when  the  first  point  is  encountered 
again.  I.e.,  L(n+1)  =  L(NMIN)  . 

Otherwise,  go  to  step  2. 

End  of  algorithm  Package  Brapping 

We  used  Sedgewick's  Pseudo  Angle  for  finding  the 
smallest  angle  in  step  3,  which  is  coded  as  the  'THETA' 
function.  This  function  returns  a  real  number  between  0-0  to 
4.0  that  is  net  the  angle  made  ny  LI  and  L2  with  the  hori¬ 
zontal  but  which  has  the  same  order  properties  as  the  true 
angle.  If  dx  and  dy  are  the  delta  x  and  y  distances  from 
some  node  to  the  anchor  node,  then  the  angle  needed  in  this 
algorithm  is  arctangent  dy/dx.  However,  the  arctangent 
function  is  likely  to  be  slow  and  it  leads  to  at  least  two 
annoying  extra  conditions  to  compute  :  whether  dx  is  zero, 
and  which  quadrant  the  point  is  in. 

In  this  algorithm  we  only  need  to  be  able  tc  compare 
angles,  not  measure  them.  Thus  it  makes  sense  to  use  a  func¬ 
tion  that  is  much  easier  to  compute  than  the  true  angle  but 
has  the  same  ordering  properties  as  the  true  angle.  A  good 
candidate  for  such  a  function  is  simply  dy  /  (dy  ♦  dx)  . 
Testing  for  exceptional  conditions  are  still  necessary,  bet 
simple  r. 

Function  THETA  (  Pseudo  Angle  ) 

Input  :  dx,uy  (delta  x  and  y  distances  from  some 
node  to  the  anchor  node) . 
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Output  i  Pseudo  angle  made  by  11  and  L 2  with  the 
horizontal  line. 

begin 

dx  =  x(L2)  -  x(L1)  :  ax  =  abs(ax)  : 
dy  =  y(L2)  -  y(Ll)  :  ay  =  abs(ay)  : 
if  (  dx=0)  and  (  dy=Q  )  then  t  =  0.0 

else  t  =  dy  /  (ax  ♦  ay  ) 
if  dx  <  0  then  t  =  2.0  -  t 

else  if  dy  <0  then  t  =  4.0  +  t 

end 

End  of  function  THETA 

Figure  2.2  shows  how  the  hull  is  discovered  in  this 
way.  He  used  Sedgewick's  pseudo  Angle  for  finding  the 
greatest  angle  selection  point. 

The  data  for  our  test  problems  is  given  in  the 
Appendix.  The  computational  results  are  summarized  in  Table 
I.  As  can  be  seen  in  Table  I,  CCAO  is  faster  than  CCCO  on 
the  small-scaled  test  problems  (below  30  nodes  ) ,  but  CCCO 
is  faster  than  CCAO  on  the  moderately  large  sized  problems 
(over  50  nodes) .  Generally,  the  accuracy  is  almost  iden¬ 
tical  in  both  cases. 

Stewart  [Ref.  6]  showed  that  in  a  large  scaled 
problem,  the  CCAO  algorithm  outperforms  any  other  insertion 
and  selection  algorithms.  Thus,  we  are  highly  motivated  to 
use  a  modification  of  the  CCAO  algorithm  for  solving  the 
time-window  constrained  TSP. 


TABLE  I 

CCHPUTATI01AL  BESDLTS  OE  CCCOr  CCAO 


CCCO  CCAO 


Prob lem 
Number 

Number 

of 

Nodes  n 

Best 

Known 

Solution 

% 

Over 

Best 

CPU 

Time 

(sec) 

X 

Over 

Best 

CPU 

Time 

(sec) 

Tu 

16 

66.603  9 

0.00 

0.0133 

0.00 

0.0066 

[2] 

22 

469.0288 

0.00 

0.0233 

0.00 

0.0100 

[3] 

22 

278.4371 

0.00 

0.0166 

0.00 

0.0066 

[5] 

51 

429.  7000 

2.  72 

0.  1897 

3.  94 

0.2562 

[6] 

76 

552.9000 

1.  64 

0.5857 

1.54 

0.6889 

III.  JjLS  ISP  WITH  HAHD  TIME  WINDOW  C0N5TBAIHTS 

A.  INTBODUCIIOH 

The  first  time-constrained  TSP  we  consider  is  the  case 
in  which  late  arrivals  are  not  allowed,  and  early  arrivals 
must  wait  for  the  opening  of  the  time  window  before  they  can 
begin  to  service  a  customer.  This  is  called  the  hard  time 
window  case  and  it  is  illustrated  in  Figure  3.1. 


The  hard  time  window  case  corresponds  to  military  opera¬ 
tions  and  to  some  civilian  distribution  problems.  Meeting  a 
deadline  rs  considered  a  critical  factor  in  this  case.  The 
soft  time  window  case  will  be  discussed  in  the  nextChapter. 

Consider  a  graph  G  =  {N,A}  composed  of  a  set  of  nodes  N 
and  a  set  of  arcs  A  connecting  these  nodes.  We  now  define 
some  notation  to  be  used  throughout  cur  discussion  of  the 
time-window-ccnstrained  TSP. 
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1  =  Lover  bound  on  the  time  window  at  node  i 

i 

(early  allowable  arrival  time  at  city  i) . 

u  =  Upper  bound  on  the  time  window  at  node  i 
i 

(latest  allowable  arrival  time  at  city  i)  . 

d  =  Time  required  to  spend  at  node  i. 
i 

(service  time  at  city  i)  . 

SPEED  =  Constant  speed  at  which  the  vehicle 
travels. 

dist  =  Distance  from  i  to  j. 
i  j 

c  =  Travel  time  from  i  to  j. 

ij 

Note  :  c  =  dist  /  SPEED. 

ij  ij 

We  use  c  and  c(i,  j)  interchangeably. 

ij 


depot  =  Depot  (home)  node. 

L  =  (  L  ( 1)  ,L(2)  ,.  L(n) 

=  A  tour  with  n  stops  visited  in  the  order 

1(1), 1(2)  ,...,L(n)  . 

ARE VT  =  Arrival  time  at  city  i. 
i 

WAIT  =  Waiting  time  at  node  i  for  the  hard  time 
i 

window. 


We  also  use  l(i),  u(i),  d(i),  ARRVT(i),  WAIT(i)  and 

1  ,  u  ,  d  ,  AERVT  ,  WAIT  interchangeably, 
i  i  i  i  i 


B.  HEURISTIC  SOLUTIOi  TECHB IQUES  FOB  BARD  TIME  WINDOWS 
1-  Nearest  Neighbor  • 

The  following  is  a  Nearest  Neighbor  heuristic 
similar  to  the  one  used  in  the  unconstrained  ISP.  At  each 
iteration  we  add  a  new  node  to  the  end  of  the  subtour.  It 


is  the  first  node  that  can  be  visited  from  the  last  node  of 


the  subtour,  taking  into  account  any  waiting  time  that  might 
be  necessary  due  to  the  lower  time  window  bounds. 

Algorithm  :  Hearest  Height or 

Input  :  Number  of  nodes,  x  and  y  co-ordinates  of  all 
nodes,  time  windows  for  all  nodes. 

Output  :  Ordered  list  of  tour,  total  travel  time. 

Step  1  .  (Initialization) 

Start  at  the  depot. 

Let  i=depot,  N'  =  {i} . 

Step  2  .  Compute  ARRVT  for  all  nodes  k€  N-N'  if  k  can 

te  visited  directly  after  i  : 

AfiRVT  =  max  (  AfiRVT  ,  1  }  ♦  d  ♦  c 

k  i  i  i  ik 

Step  3  .  If  AfifiVT  >  u  ,  then  stop  ('  no  feasible 

k  k 

solution' ) . 

Step  4  .  If  AfifiVT  <  1  ,  then  cost  =  1 

k  k  k  k 

Otherwise,  cost  =  ARRVT 

k  k 

Step  5  .  (  Nearest  Neighbor  Selection) 

Choose  the  node  k€  N-N'  such  that  cost  is 

k 

a  minimum.  I.e,  find 

k  =  argmin  cost  s.t.  j  6N-N'. 

3 

Step  6  .  (Insertion) 

Insert  k  after  i,  add  k  to  subtour  N ' ,  and 
let  i  =  k. 

Step  7  .  If  N '=  N ,  go  to  next  step. 

Otherwise,  return  to  step  2. 

Step  8  .  Compute  total  travel  time,  then  stop. 

Total  travel  time  =  max  {  AfifiVT  ,  1  } 

k  k 

♦  d  ♦  c 

k  k,  depot 

End  of  algorithm  Nearest  Neighbor 


Tiiis  solution  was  constructed  by  starting  at  the 
depot  and  moving  to  the  nearest  neighboring  customer  that 
has  not  yet  teen  visited  as  long  as  the  upper  bound  level 
was  not  violated.  This  heuristic  may  fail  to  solve  the 
problem . 

2.  SCCO 

This  algorithm  is  designed  for  the  case  when  some  of 
the  nodes  do  not  have  time  windows.  He  call  these  nodes  *' 
time  free  n . 

SCCO  is  similar  to  the  cheapest  selection,  cheapest 
insertion  method  for  the  unconstrained  TSP,  except  that  the 
nodes  with  time  windows  are  treated  differently  from  the 
time  free  nodes.  The  nodes  with  windows  are  inserted  in 
order  of  increasing  upper  time  window  bound. 

The  time  free  nodes  are  inserted  between  these  nodes 
by  cheapest  selection  and  cheapest  insertion,  for  as  long  as 
the  upper  bound  time  window  will  allow.  In  the  end,  a 
Modified  oropt  is  used  to  improve  the  solution. 

There  is  one  possible  difficulty  with  this  approach. 
It  may  become  impossible  to  reach  some  of  the  time- 
constrained  nodes  before  their  upper  bound.  In  this  case, 
we  must  delete  some  node  (s)  from  the  subtour.  Whenever  we 
see  an  upper  oouna  that  cannot  be  satisfied,  we  select  a 
node  to  delete  by  the  following  criteria. 

The  first  criterion  is  the  width  of  the  time  window. 
Hence,  time-free  nodes  are  considered  first.  Then,  if 
several  nodes  in  the  subtour  are  tied  for  the  widest  time 
window,  we  select  for  deletion  the  node  tnat  results  in  the 
greatest  tine  saved.  The  algorithm  is  summarized  as 
fellows . 

Algorithm  :  Successive  Cheapest  Cheapest  Oropt  (SCCO) 


Input  :  cumber  of  nodes,  x  and  y  co-ordinates  of  all 


Output  : 


Step  1  . 


S  tep  2  . 


nodes,  time  windows  for  all  nodes, 
ordered  list  of  tour,  total  travel  time- 

(Initialization) 

Start  at  the  depot, 
let  i=depot,  N*  =  {i}. 

Set  k  =  argmin  u  s-t.  j  €  N-N*. 

j 

If  k  is  time  free  node,  then  set  k  =  depot. 


Step  3  .  Calculate  AERVT  . 

k 

A  HR  VT  =  max  {  ARRVT  ,  1  }  ♦  d  *■  c 

k  i  i  i  ik 

Step  4  .  If  ARRVT  <  u  ,  then  go  to  step  5. 

k  k 

Otherwise,  select  time  free  node  m  €  N*  which 
results  in  the  greatest  time  saved  for  dele¬ 
tion.  Delete  node  m  from  N ' ,  go  to  step  3. 
Step  5  .  Add  node  k  to  the  subtour  N  • . 

Step  6  .  Insert  time  free  node  j£  N-N'  between  nodes  i 
and  k  by  cheapest  insertion  and  cheapest 
selection  (same  as  CCCO)  until  ARRVT^does  not 

exceed  u 


If  ARRVT  <  1  ,  then  set  ARRVT  =  1  . 

k  k  k  k 

Step  7  .  If  M'=N,  then  go  to  next  step. 

Otherwise,  let  i  =  k  and  go  to  step  2. 

Step  8  .  Apply  the  Modified  Oropt  procedure  to  the 
current  tour.  Stop  when  no  further  improve¬ 
ments  can  be  found. 

End  of  algorithm  SCCO 

"Successive"  means  select  the  node  successively  by 
the  smallest  upper  bound.  In  the  SCCO  algorithm,  if  the 


salesman  arrives  before  the  lover  bound  of  the  time  window, 
adding  waiting  time,  we  set  the  arrival  time  equal  to  the 
lover  bound. 

The  Modified  Oropt  procedure  for  improving  the  solu¬ 
tion  is  described  below.  This  procedure  consider  only  those 
exchanges  that  would  result  in  a  node  being  inserted  between 
two  other  nodes  in  the  current  tour. 


Figure  3.2  Current  Tour  before  Modified-Oropt. 

Figures  3.2  and  3.3  are  helpful  to  understand  how 
the  procedure  works.  In  both  figures,  i,j,k,l,  and  m  are  the 
nodes  in  the  current  tour.  Nodes  1  and  m  are  considered  to 
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be  adjacent  tc  node  k.  A  test  is  then  conducted  to  deter 
mine  if  node  k  can  be  located  between  two  other  nodes,  sue 
as  i  and  j,  sc  that  it  results  in  reduced  total  travel  time 
If  it  can,  we  make  the  appropriate  arc  exchanges,  the 
update  the  total  cost  and  route  orders. 


Figure  3.3  Improved  Tour  after  Modified-Oropt. 


In  this  example,  the  three  arcs  (i, j) ,  fk,l) ,  and 
(k,m)  are  removed  and  replaced  by  (i,_k),  (j,k),  and  (l,m) 

When  no  further  exchanges  improve  the  solution,  the  algo 
rithm  terminates. 


Figure  3.4  Subtour  in  SCCO  Procedure 


Figures  3. 4  -  3.6  illustrate  the  SCCO  algorithm  for 
the  TSP  with  hard  time  windows  given  in  Appendix  F  as  in 
test  problem  [1  ].  In  this  problem  10  of  16  nodes  have  time 
windows.  The  ether  6  nodes  are  time  free. 

First,  the  subtour  starts  at  the  depot  (node  16)  and 
we  insert  the  node  with  the  smallest  upper  bound  (node  12). 
He  examine  all  nodes  which  could  be  inserted  between  16  and 
12  as  long  as  the  upper  bound  on  node  12  is  observed.  In 
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Figure  3.5  Intermediate  Subtour  in  SCCO  Procedure. 


this  case  there  is  no  such  node.  Then  we  select  the  next 
smallest  upper  bound  (node  1 4) ,  add  it  to  the  tour,  look  for 
nodes  to  insert  before  it,  and  continue  in  this  manner.  Now 
we  have  formed  the  partial  tour  16,  12,  14,  11,  6,  3  as 
shown  in  Figure  3.4. 

As  shewn  in  Figure  3.5,  we  can  insert  3  time  free 
nodes  between  node  6  and  node  3.  These  insertions  are  made 
according  to  the  cheapest  insertion  and  cheapest  selection 


Figure  3.6  Final  Subtour  for  SCCO 


criteria.  He  do  not  make  any  further  insertions  because 
they  would  cause  a  time  window  violation  at  node  3. 

Figure  3.6  shows  the  final  tour  for  the  SCCO 
heuristic.  This  tour  is  passed  to  a  Modified  Oropt,  but  in 
this  case  it  will  find  no  improvement. 


3.  SCAO 


This  heuristic  is  identical  tc  SCCO  except  for  the 
use  of  the  greatest  angle  selection  criterion  for  the  time- 
free  nodes,  instead  of  cheapest  selection. 


Algorithm  :  Successive  Cheapest  Angle  Oropt  (SCAO) 


Input  ;  Number  of  nodes,  x  and  y  co-ordinates  of  all 
nodes,  time  windows  for  all  nodes. 

Output  :  Ordered  list  of  tour,  total  travel  time- 


step  1 


S  tep  2 

Step  3 


S  tep  4 


(Initializ  ati on) 

Start  at  the  depot. 

Let  i=depot,  N*  =  £i } . 

Set  k  =  argmin  u  s.t.  j  €  N-N’. 

j 

If  k  is  time  free  node,  set  k  =  depot. 

Calculate  ARE  VT  . 

k 

AEfiVT  =  max  {  ARRVT  ,  1  }  +  d  +  c 

k  i  i  i  ik 

If  ARRVT  <  u  ,  then  go  to  step  5. 
k  k 


S  tep  5 
S  tep  6 


Otherwise,  select  time  free  node  m€N'  which 
results  in  the  greatest  time  saved  for  dele¬ 
tion.  Delete  node  m  from  N',  go  to  step  3. 

Add  node  k  to  the  subtour  N'. 

Insert  time  free  node  j  €  N-N'  between  nodes  i 
and  k  by  cheapest  insertion  and  greatest 
angle  selection  (same  as  CCAO)  until  ARRVT 

k 


does  not  exceed  u 

k 

If  ARRVT  <  1  ,  then  set  ARRVT  =  1  . 

k  k  k  k 


Step  7  .  Let  i  =  k- 

If  N'  =  N,  then  go  to  next  step. 
Otherwise,  go  to  step  2. 


Step  8  .  Apply  the  Modified  Oropt  procedure  to  the 
current  tour.  Stop  when  no  further  improve¬ 
ment  can  be  fcund. 

End  of  algorithm  SCAO 

This  algorithm  is  same  as  SCCO  except  greatest  angle 
selection  is  used  instead  of  cheapest  selection,  as  in  SCCO. 

4.  SLACJ5 

This  heuristic  was  suggested  by  Professor  Rosenthal. 
It  is  designed  for  the  case  when  the  widths  between  the 
upper  and  lower  bounds  of  the  time  windows  are  relatively 
large. 

In  this  heuristic,  the  SLACK  is  the  most  important 
concept.  The  SLACK (i)  can  be  defined  as  the  maximum  amount 
of  time  by  which  arrival  at  node  i  can  be  delayed 
without  causing  an  upper  bound  to  ne  violated  for  a  node 
currently  on  the  tour. 

The  SLACK  function  can  be  defined  as  a  recursive 
function  as  fcllcws. 

SLAC  K  (L  ( i)  )  =  min  {  u(L(i))  "  ARBV  T  (L  (i) )  , 

SLACK  (L  (i+  1) )  +  W  AIT  (L  (i) )  } 

where 

HAI T  (L  (i  )  )  =  max  {0,  1(1  (i))  -  AREVT  f  L  (i )  )  } 

The  first  element  of  this  recursive  function  is  the 
difference  between  the  upper  bound  and  arrival  time  at  node 
l(i).  The  second  one  is  the  sum  of  next  node's  SLACK  and 
waiting  time  cf  node  L(i).  The  minimum  of  these  two 
elements  is  a  possible  delay  time  of  the  arrival  time  at 
ncde  L(i)  without  violating  the  upper  bound  of  all  nodes 
after  L  (i)  in  the  current  tour. 


The  advantage  of  this  recursive  function  is  that  it 
is  easy  to  calculate  a  possible  delay  time  without  calcu¬ 
lating  new  arrival  times  for  all  nodes  after  L  (i) .  The  algo¬ 
rithm  is  summarized  as  follows. 

Algorithm  :  SLACK 

Input  :  Number  of  nodes,  x  and  y  co-ordinates  of  all 
nodes,  time  windows  for  all  nodes. 

Output  :  Ordered  list  of  tour,  total  travel  time. 

Step  1  .  (Initialization) 

Start  at  the  depot.  Let  N'  =  {depot}. 

Step  2  .  Sort  the  upper  time  windows. 

U  =  (  u  ,u  , - - u  ) 

1  2  n 

s.t.  u  < u  < . . .< u 
1  2  n 

Step  3  .  Set  k  =  argmin  u  s.t.  i£  N-N*. 

i 

Step  4  .  Find  a  node  L(ISTAR)  after  which  node  k 
can  be  inserted  in  the  current  sequence, 
if  such  a  node  exists.  Go  to  step  7. 

(The  criteria  by  which  we  determine  xf  an 
insertion  can  be  made  are  given  below.) 

Step  5  .  If  there  is  no  such  place  to  insert  node  k, 
then  try  to  find  a  node  L(ISWAP)  in  the 
current  sequence  such  that  k  can  replace 
L(ISWAP)  and  I  (IS HAP)  has  a  good  chance  of 
being  reinserted  somewnere  else. 

Select  ISWAP  which  has  the  largest  time 
window  width  among  candidates  for  ISWAP. 

If  there  is  no  candi-date,  then  stop. 

(  '  no  feasible  solution  '  ) 

.  Do  swap  (  add  k  to  N',  and  delete  L  (ISWAP) 
from  N',  and  set  k  =  L(ISWAP)  ), 
then  update  slack  and  arrival  times. 


S  tep  6 


Step  7  .  Select  the  node  which  results  in  the  minimum 
additional  travel,  i  .e,  the  node  k  which 
minimizes  the  following  guantity. 
c  (L  (I)  ,  k)  ♦  c(k,L(I+1))  -  c(L(I)  ,L<1+1))  . 

Step  8  .  Insert  k  after  L(ISTAfi)  ,  and  add  k  to  N', 
and  update  slack  and  arrival  times. 

Step  9  .  If  N'=N,  then  stop. 

Otherwise,  go  to  step  3. 

End  of  algorithm  SLACK 


This 

pr ocedure 

starts 

with 

sorting 

an  array 

", ' 

u  ,  .  .  .  , 

2 

u  into 

n 

ascending 

order 

using  a 

heapsort 

[  Ref. 

20], 

This  u 

array  is 

used  to 

select  a 

node  k  in 

ascending  order  for  insertion.  Since  the  upper  bound  cannot 
be  violated,  this  step  is  performed.  Then  find  a  node  1(1) 
after  which  ncde  k  can  be  inserted  in  the  current  sequence, 
if  a  such  a  ncde  exists. 

There  are  two  tests  which  must  be  administered  to 
determine  if  k  can  be  inserted  after  L(I).  First,  the 
arrival  time  at  node  k  if  k  succeeds  L(I),  which  is  called 
TEST1  must  not  be  greater  than  the  upper  bound  u  .  Second, 
if  k  precedes  1(1+1),  tnen  the  resulting  delay  in  arrival  at 
1(1+1),  which  is  called  TEST2  ,  must  not  greater  than 
SLACK  (1  (1+1 ))  .  We  can  calculate  TESJ1,  TEST2  as  follows. 

TEST1  =  Arrival  time  at  node  k  if  k  succeed  i(i). 

=  max  {AfiRVI  (L  (i)  ) ,  1  <1  (i) )  )  +  d(l(i))  +  c(l(i),k). 

TEST2  =  Delay  in  arrival  at  L(i*1)  if  k  precedes  L(i+1). 

=  max  {  IEST1 ,  l(k)  }  +  d  (k )  ♦  c(k,L(i+1)). 


If  there  exists  more  than  one  node  L(I)  after  which 
node  k  can  be  inserted,  we  select  L  (I)  according  to  the 
criterion  of  least  additional  travel  time.  This  additional 
travel  time,  called  TEST3,  is  given  by 

TEST3  =  c(l|I),k)  ♦  c{k,  1(1*1))  +  C(1(I)  ,L(I+1) )  . 

When  we  insert  node  k  after  L(I)  ,  we  update  the 
arrival  times  and  SLACKS.  In  the  updating  process,  we 
compute  updated  values  of  SLACK  only  for  the  nodes  whose 
SLACK  actually  changes  as  a  result  of  the  insertion. 

If  there  is  no  place  to  insert  node  k,  we  call  a 
subroutine  called  'TSWAP'.  TSKAP  tries  to  find  a  node 
L  (ISWAP)  in  the  current  seguence  such  that  k  can  replace 
L  (I SWAP )  and  L(IS»iAP)  has  a  good  chance  of  being  reinserted 
somewhere  else.  TSWAP  uses  TEST1  ans  1EST2  to  find  a  candi¬ 
date  for  IS  HAP  and  then  uses  a  largest  time  window  width  to 
select  ISWAP.  If  there  exists  such  a  ISWAP,  then  we  do  the 
swap  and  update  SLACKS  and  arrival  times  and  try  to  insert 
again. 

C.  EXACT  SOLUTION  TECHIIQUES  FOB  BABD  TINE  WINDOWS 

1.  State-Space  Reiaxati on  Procedure 

A  dynamic  programming  model  of  the  time-constrained 
TSP  has  been  developed  by  Christofides  et  al.  [Ref.  5].  we 
applied  their  approach  to  compute  bounding  information 
within  a  branch  and  bound  algorithm. 

Consider  the  TSP  defined  on  the  graph  G  =  (N,A)  with 
the  time  window  constraints,  where  N  is  a  set  of  all  nodes 
of  G,  and  A  is  a  set  of  arcs.  Let  R(j)  be  the  set  of  nodes 
from  which  it  is  possible  to  go  directly  to  node  j.  We  can 
initially  set  R(j)  =  N  -  [i  |1;  ♦  dj.  ♦  c-j  >  uj  } ,  because 
it  is  impossible  to  go  directly  from  ncde  i  to  node  j  if  the 
earlist  possible  arrival  time  at  node  j  exceeds  the  upper 
time  window  cf  node  j. 


46 


Let  f(S,j)  be  the  duration  of  the  least  time  path 
starting  at  node  1  passing  through  every  node  of  SSS'=N-{1} 
and  finishing  at  node  j.  For  a  given  S  and  j,  xe  can  calcu¬ 
late  a  minimum  arrival  time  in  node  j  as 


T(S#j)  =.  min  r  f(S-j,i)  ♦  d .  ♦  c.J. 

ie(s-gnaa))  i  13 


Then, 


f(S,j)  = 

T(S,j)  , 

if 

1  .  <  T  (S  ,  j ) 

- 

1 

3 

00  , 

if 

TlS,j)  <  1 

2 

if 

T  (S,  j)  >  u 

the  initialization: 

f  ({j},  3) 

=  C  (1,  j) 

9 

if  1 ,  <  c 
if  c"5  <  1 

=  1 

9 

j 

1  j 

—  00 

9 

if  c  >  u 

1j 

(3.  1) 


In  equation  (3.1)  the  minimum  arrival  time  in  node  j  passing 
through  the  nodes  in  the  set  S  can  be  described  as  the  sum 
of  three  terms  ;  the  first  is  the  duration  of  the  least  time 
path  passing  through  the  nodes  in  the  set  S- {j}  and  ending 
in  node  i,  the  second  is  the  time  required  to  spend  in  node 
i,  and  the  third  is  the  travel  time  from  node  i  to  node  j. 

The  f(S,j)  can  be  calculated  for  all  subsets  5  of  S' 
and  for  all  nodes  j  by  using  eguation  (3.1)  recursively. 
Finally,  the  optimum  solution  can  be  calculated  as 

min  [  f(S',i)  ♦  d  *■  c  ]. 
i  e  S'  i  il 

Since  the  computer  storage  requirements  increase 

exponentially  with  the  size  of  the  problem,  this  method  is 

limited  to  small  problems.  The  total  number  of  f  (S,  j)  ,  when 

/n-l  \ 

S  contains  k  nodes,  is  k(  v  1,  since  f(S,j)  must  be  calcu¬ 


lated  for  all  subsets  S  o 


and  since  each  node  in  S  must 


be  considered  as  a  possible  end-node  3.  Therefore  the 
storage  requirement  for  f  (S, j)  in  a  n  node  problem,  is  given 
by  [Bef  .  21  ]. 


n-1  /  n-1 \  n-: 

m  k  ( * )  ■  2  • 


(3.2) 


The  storage  requirements  to  solve  a  22  node  problem 
exceed  22,020,096.  For  relaxing  this  limitation, 
Christofides  et  al.  [Ref.  5]  proposed  a  state  space  relax¬ 
ation  procedure  which  is  analogous  to  Lagrangean  relaxation 
[Bef.  22]  in  integer  programming.  The  state  space  associ¬ 
ated  with  a  given  dynamic  programming  recursion  is  relaxed 
in  such  a  way  that  the  solution  to  the  relaxed  recursion 
provides  a  bound  which  could  be  embedded  in  a  general  branch 
and  bound  method  [Bef.  23].  He  describe  Christofides  et  al*  s 
method  for  doing  this  below. 

Consider  the  dynamic  programming  formulation  (3.1) 
The  state  variable  in  that  formulation  is  (S,  j)  ,  and  the 
stage  is  the  cardinality  of  S.  Let  g  (S)  be  a  mapping  from 
the  domain  of  (S,j)  to  some  other  vector  space  (y(S),j). 
Let : 


h  (g(S)  ,  j)  =  {  (g  (s-j)  ,i)  |  i  e  (s-jn  8(3))  } 


(3.3) 


Since  we  are  interested  in  lower  bounds  to  the  TSP 
with  time  constraints,  H(g(S),j)  in  (3.3)  may  be  replaced  by 
any  larger  set  that  is  easier  to  compute.  Thus,  H(g(S),j) 
can  be  defined  by  the  following  eguaticn: 


H  (g  (S)  ,  j)  =  {  (g  (S-j)  ,  i)  |  i  e  E  ( g  (S) ,  j)  } 
where  (S-j  n  fi(  j))  i  E  (g(S) ,  j) . 


(3.4) 


Fcr  calculating  the  lower  bound  of- the  problem,  equation 
(3.1)  can  be  changed  to  the  following  equation: 

T  <g(S),  j)  =  min  [f  (g(S-j)  ,i)  *  d.  ♦  c.  .]  (3.5) 

(g(S-g) ,i)  e  H ( g ( Sf ,1  i  il 


nn 


i  e  E  (g(S),  j) 


[f  (g(S-j)  ,i)  ♦  d,  ♦ 


c.  .3 

13 


(3.6) 


This  gives  us: 

f(g(S),j)  =  T(g(S),j),  if  l.<  T(g(S),j)  <u. 

=  1  .  ,  if  T?g(S)  ,j)  <  1  .  3 

=  «?  *  if  T  (g  (S)  ,  j)  >  u3. 

3 

With  the  initialization: 


f  <g(j)  #  3)  =  c 
=  1 


,  if  1  <  c  <  u 

lj  j  .1ji  j 

,  if  c  <1 

3  .  ,  lj  3 

,  if  c  >  u 
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Finally,  tne  optimum  solution  can  be  calculated  as 

min  [  f(g(S')/i)  ♦  <1  ♦  c  ]. 

i  e£(g(()  ,1)  i  i  1 


The  mapping  can  be  selected  frcm  any  separable  func¬ 
tion.  Chnstofides  et  al.  used  the  following  mapping 
function. 

g  (S)  =  |  S|.  (3.7) 

Then  eguation  (3.6)  becomes  : 

T(k,j)  =  min  [  f(  k-1,i  )  ♦  d  ♦  c  ]  (3.8) 

i  e  E (k, j)  x  i j 

where  k  =  | S  J  >  1 . 


This  gives  us: 

f  (k,j)  =  Kk,j)  ,  if  1.5  T (k, j)  $  u  , 

=  1  .  r  if  T^krj)  <  1  .  3 

=  co3  ,  ifT(k,j)>u3. 

3 


With  the  initialization: 


f  Mr  j) 


C (1, j) 


9 


if  1  <  c 


<  u 


Finally*  the  optimum  solution  can  be  calculated  as 

.  Bin  [  f(lS'|,i)  ♦  d.  ♦  c .  ]. 

1  e  E  ( |  N  | ,  1)  i  il 


2.  Add  it ional  Condition 

In  the  previous  section,  we  discussed  Christofides 
et  al. * s  state  space  relaxation  procedure  which  provides  a 
lower  bound  on  the  TSP  by  reducing  a  state  space  in  dynamic 
programming.  This  lower  bound  is  effective  in  branch  and 
bound  only  if  it  is  a  tight  bound.  This  is  similar  to  the 
case  in  integer  programming  where  the  effectiveness  of 
Lagrangean  relaxation  in  producing  bounds  is  relative  to  the 
integer  programming  formulation.  A  redundant  state-space 
condition  can  be  helpful  to  get  a  better  bound-  For  this 
purpose,  an  additional  condition  was  used  by  Christofides  et 
al.  to  avoid  loops  formed  by  three  consecutive  nodes  [Ref. 
5].  This  can  be  done  in  the  following  way. 

Let  k  =  1 S | .  Let  f  (k,j,1)  be  the  duration  of  the 
least  time  path  from  tne  initial  state  to  state  (k,j) 
without  loops  formed  by  three  consecutive  nodes.  Let 
f(k,j,2)  be  the  duration  of  the  second  least  time  path  from 
the  initial  state  to  state  (k,j)  without  loops  formed  by 
three  consecutive  nodes.  Let  p(k,  j,m)  be  the  predecessor  of 
j  on  the  path  corresponding  to  f(k,j,m).  With  the  above 
definition,  recursion  (3.8)  becomes: 


T(k,  j,  1)  =  min  [  f(K-1,i,m)  ♦  d  ♦  c  1, 

l  e  E(k,  j)  i  ij 


(3.9) 


where  m  =  1,  if  p<k-1,i,  1)*j 
=  2,  otherwise. 

This  gives  us: 


f(k,j,1)  =  I(k,j,1),  if  1  <  T  (k,  j,  1)  <  u 

=  1.  ,  if  T"j k,  j ,  1)  <  1.  1 


(3.  10) 


00 


if  T(k,j,1)  >  u 

3 


t 


Hecursion  for  f  (k,  j,2)  can  be  written  in  the 
following  way. 

(3.11) 


(3-12) 


(3.13) 

(3.14, 
(3.15) 

Since  the  additional  condition  can  avoid  considera¬ 
tion  of  a  useful  lower  bound,  we  considered  f(k-1,i,2)  in 
recursion  (3.9,  and  (3.11)  only  when  the  predecessor  of  i  on 
the  path  corresponding  to  f(k-1,i,1)  is  j.  If  we  do  not 
consider  the  second  least  time  path  in  case  of  p  (k-  1,  i,  1)  =  j, 
then  f(q(5),j)  dees  not  guarantee  the  lower  bound  of  f(S,j). 

For  this  example,  let's  consider  a  4  node  TSP  with 
time  constraints.  Node  A  is  the  starting  node.  D  is  the 
time  free  node.  The  lower  bound  of  node  B  is  9,  the  upper 
bound  of  node  B  is  11,  the  lower  bound  of  node  C  is  19,  and 


let : 


T(k,j,2)  =  .  min 


min  f  f(k-1,i,m)  +  d  +  c  ], 

.  ieEU^jl  i  ij 

l  *p(k,3,1) 


where  m  =  1,  if  p(k-1,i,1)*j 
=  2,  otherwise. 

Ibis  gives  us: 

f  (k,  j,  2)  =  T(k,  j,2)  , 

=  1 


if  1  <  T  (k,  j,  2)  <  u 

if  T^k,  j,  2)  <1.  J 

if  T (k, j, 2)  >  u3. 

3 


The  initialization  is 


f(1,i,1)  =  c ( 1 , i,  ,  if  1.  < 

i 

=  1 

i 

=  00 


*  if 

#  if 


c  <  u 
1i  i 
c  <1 

1i  i 

c  >  u 

1i  i 


and 


f  ( 1#i»  2)  =  « 

Finally,  the  optimum  solution  can  be  calculated  as 

min  [  f(|S'|,i,1)  ♦  d  +  c  ]. 
i  e  E ( | S | ,  1 )  i  i  1 
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tie  upper  bound  of  node  C  is  21.  Suppose  service  time  at 
each  node  is  zero.  Figure  3.7  shows  an  optimal  route  for 
this  prcbiem. 


From  equation  (3.13)  we  can  get: 
f  (1,B,  1)  =  10, 

f(1,C,1)  =  19, 

f  (1,0,  1)  =  7.07 

Now  applying  equation  (3.9)  recursively  with  i=  1 ,  for  k  =  2  we 
can  get: 
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f  <2,B,1) 


oo  r 


f(2,C,1)  =  19, 
f  (2  ,  D  ,  1  )  =  17.07 

Similarly,  for  k=3 

f  (3,  B  ,  1 )  =  oo  , 
f  (3,C,  1)  =  =o  , 
f  (3,0,  1)  =  “  - 

We  can  see  easily  that  f(3,D,1)  is  not  a  lower  bound  of 
f  ({3.C,£}  ,D )  . 

3-  Branch  and  Bound  Procedure 

In  this  section  we  introduce  branch  and  bound 
enumeration  which  is  used  to  eliminate  subtours  in  the  solu¬ 
tion  of  the  state  space  relaxation  procedure.  Since  the 
state  space  relaxation  procedure  is  a  relaxation  of  the  TSP 
with  time  constraints,  the  solution  to  the  state  space 
relaxation  procedure  provides  a  lower  bound  on  the  optimal 
value  of  the  1SE  with  time  constraints.  Any  heuristic  solu¬ 
tion  can  provide  an  upper  bound.  He  denote  some  notation  to 
explain  this  algorithm  as  follows. 

FLBD  =  The  lower  bound,  which  is  the  optimal  solut¬ 
ion  to  the  state  space  relaxation  procedure, 
on  the  optimal  solution  to  the  TSP  with 
time  constraints  given  restrictions  at  the 
current  node. 

Z  =  Current  upper  bound. 

STACK  =  Array  which  represent  decision  tree.  It  con¬ 
tains  arc  lists  which  have  the  same  head  in 
optimal  tour  to  the  state  space  relaxation 
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procedure  given  restrictions  at  the  current 
node. 

[ c»  ]  =  Travel  time  matrix  given  restrictions  at 

ij 

the  current  node. 

There  are  two  types  cf  tree  search.  One  is  depth- 
first  search,  the  other  is  breadth-first  search  [Eef.  24]. 
lie  used  depth-first  search  since  breadth-first  search 
required  substantially  more  storage.  Depth-first  search 
simply  means  that  when  a  separation  is  defined,  one  of  the 
nodes  created  by  the  separation  is  immediately  selected  to 
be  the  next  suiproblem,  and  when  a  node  is  fathomed,  the 
enumeration  always  backtracks  to  the  most  recently  created 
live  node- 

One  of  the  most  important  requirements  of  any  branch 
and  bound  algorithm  is  tight  bounds.  The  closer  the  bounds 
are  to  the  optimal  solution,  the  fewer  nodes  must  be  enumer¬ 
ated.  We  used  the  SCCO  heuristic,  which  was  described  in 
section  B.2,  as  an  initial  upper  bound.  The  lower  bound  is 
obtained  from  eguation  (3.15). 

To  save  computing  time  we  need  a  criterion  to  decide 
whether  or  not  the  branching  should  be  continued.  If  FLBD 
is  greater  than  Z,  then  the  node  is  fathomed  since  explicit 
enumeration  need  not  be  extended  belcw  the  current  node. 
For  branching  we  consider  the  arcs  which  have  the  same  head 

node  in  the  directed  graph  since  each  arc  must  have  a 

different  head  in  the  TSP  solution.  If  there  is  no  such 
arc,  then  that  solution  is  a  feasible  solution.  After  all 
nodes  of  the  tree  are  fathomed,  a  feasible  solution  which 

has  the  same  value  as  the  upper  bound  is  an  exact  solution 

to  the  TSP  with  tame  window  constraints. 

The  following  branch  and  bound  algorithm  is  used  in 
the  programs  written  for  exact  solution. 


Algorithm  :  Branch  and  Bound  Procedure 

Input  :  Iotal  travel  time  of  heuristic,  travel  time. 
Output:  Ordered  list  of  tour,  total  travel  time. 


Step  1 

Step  2 
Step  3 

Step  4 

Step  5 
Step  6 


Step  7 


(Initialization) 

Let  Z  =  the  optimal  solution  of  SCCO. 

STACK  =  empty. 

[  c\  ]  =  [  c.  .  ] 
in  13 

Compute  FLBD  given  restrictions  defined  by 

[  c'  ].  If  FLED  >  Z,  go  to  step  5. 
i  j 

(Construct  the  tree) 

Put  all  arc(i,j)  which  have  the  same  head  j 
in  directed  graph  on  STACK. 

If  there  is  no  such  arc,  save  feasible  route 
and  update  Z  =  FLBD  then  go  to  step  5. 

Let  travel  time  of  arc(i,j)  which  is  in  the 
top  of  STACK  be  infinite,  then  go  to  step  2. 
(i.e.  ,  c  '  =oo  .) 

(Backtrack) 

If  STACK  =  empty,  go  tc  step  7. 

If  travel  time  of  arc(i,j)  which  is  in  the 

top  of  STACK  is  finite,  let  travel  time  of 

that  arc  (i,  j)  be  infinite,  then  go  to 

step  2.  (i.e.,  cf  =  . ) 

ij 

Otherwise,  let  travel  time  of  that  arc(i,j) 
be  original  travel  time  of  that  arc(i,j)  and 
remove  tnat  arc(i,j)  from  top  of  the  STACK. 
Go  to  step  5.  (i.e.  ,  c*  =  c  ) 

ij  ij 

(termination) 

If  there  is  a  feasible  route,  then  the 
optimal  travel  time  =  Z. 

Otherwise,  there  is  no  feasible  solution. 


End  of  algorithm  Branch  and  Bound  Procedure 


17.  IgE  ISP  WITH  SOFT  TIME  WINDOW  CONST SAINTS 


A.  INTBODDCTIOB 

The  second  time-cons train  ted  TSP  ve  consider  is  the  case 
in  which  both  late  and  early  arrivals  are  allowed  by  paying 
a  penalty  cost.  The  penalties  are  allowed  to  be  different 
for  early  and  late  arrivals.  Tne  penalty  cost  is  calculated 
as  follcws. 

jpper  penalty  cost  =  sax  [  G,  upper  penalty  constant 

x  (  arrival  time  -  upper  bound  )  ]. 

Lower  penalty  cost  =  max  [  0,  lower  penalty  constant 

x  {  lower  bound  -  arrival  time  )  ]. 

In  fact,  the  upper  penalty  constant  is  greater  than  the 
lower  penalty  constant  in  most  cases.  Figure  4.1  may  be 
helpful  to  understand  this  case. 

This  approach  makes  every  problem  feasible,  no  matter 
what  the  time  windows  are,  i.e,  even  if  it  is  infeasible  in 
the  hard  time  window  case.  This  reflects  a  practical  point 
of  view,  especially  when  it  is  possible  to  save  a  great  deal 
of  mileage  by  allowing  a  small  amount  of  time  window 
violation. 

In  this  Chapter,  we  considered  one  unit  of  cost  to  be 
the  same  as  one  unit  of  time.  In  real  world  problems,  it  is 
possible  to  get  a  cost  by  multiplying  traveling  time  by  some 
constan  t. 

We  use  the  notation  Ip  and  up  for  the  lower  and 

k  k 

upper  penalty  cost  at  node  k. 
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lower 

penalty 

cost 


upper — > 

penalty 

cost 


early 

arrival 


lower 

bound 


upper 

bound 


time  window 
for  city  i 


— x - 

late 

arrival 


Figure  4.1  Diagram  for  Soft  Time  Window  Case. 


B.  HEURISTIC  SOlOTIOl  TECHNIQUES  FOB  SOFT  TIME  I IN DO IS 

I^afest  Neighbor 

This  heuristic  is  similar  to  the  hard  time  windows 
except  it  takes  into  account  any  penalty  cost  that  might  be 
necessary. 

Algorithm  ;  Nearest  Neighbor 

Input  :  Number  of  nodes,  x  and  y  co-ordinates  of  all 
nodes,  time  windows  for  all  nodes. 

Output  :  Ordered  list  cf  tour,  total  cost. 

Step  1  .  (Initialization) 

Start  at  the  depot. 

Let  i=depot,  N*  =  (i)  ,  cost  =  0. 

i 

Step  2  .  Compute  ARRVT  for  all-ncdes  k  €  N- N ' 

ARRVT  =  AESVT  ♦  d  -*•  c  . 
k  i  i  ik 

cost  =  cost  *  d  ♦  c 
k  i  i  ik 


58 


S  tep  3 


step  4 


If 

AREVT 

< 

1 

,  then 

cost  = 

cost  ♦ 

iP 

k 

k 

k 

k 

If 

Afifi  VT 

> 

u 

,  then 

cost  = 

cost  + 

up 

k 

k 

k 

k 

(  Nearest  Neighbor  Selection  ) 

Select  the  node  k  €  N-N '  such  that  cost 
is  a  minimum.  I.e.  ,  find 
k  =  argmin  cost  s.t.  j  €  N-N • . 

j 


Step  5  .  (  Insertion  ) 

Insert  k  after  i,  add  k  to  subtour  N'#  and 
let  i  =  k  . 

Step  6  .  If  N'  =  N  ,  then  go  to  next  step. 

Otherwise,  go  to  step  2. 

Step  7  .  Compute  total  cost,  then  stop. 

Total  cost  =  ccst  ♦  d  ♦  c 

k  k  k, depot 


End  of  algorithm  Nearest  Neighbor 

This  solution  was  constructed  by  starting  at  the 
depot  and  moving  to  the  nearest  neighboring  customer  that 
has  not  yet  been  visited.  The  term  "nearest"  is  modified  in 
the  sense  that  we  add  a  penalty  cost  to  the  travel  time  if 
the  time  window  for  city  i  is  violated. 


2.  SCCO 


This  algorithm  is  also  designed  for  the  case  when 
there  is  a  combination  of  tight  time  window  nodes  and  time 
free  nodes.  The  strict  observance  of  the  upper  bound  in  the 
hard  time  windows  is  replaced  by  a  penalty  cost. 

Algorithm  :  SCCO 

Input  :  Number  of  nodes,  x  and  y  co-ordinates  of  ail 
nodes,  time  windows  for  all  nodes. 

Output  :  Ordered  list  of  tour,  total  cost. 


Step  1  .  (Initialization) 
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S  tep  2 . 


Start  at  the  depot. 

Let  i=depot,  N •  =  {i}  ,  cost  =  0. 

i 

Set  k  =  argmin  u  s.t.  j  €  N-N • . 

If  k  is  time  free  node,  then  set  k  =  depot. 


Step  3  .  Insert  node  k  in  the  sabtour  N'. 

Compute  ARBVT 

k 

ARRVT  =  ABRVT  ♦  d  ♦  C 

k  i  i  ik 

Step  4  .  Insert  time  free  node  j£  hi— N  •  between  nodes 

i  and  k  by  cheapest  insertion  and  cheapest 

selection  {  same  as  CCCC)  until  AERVT  does 

k 

not  exceed  u  . 

k 

Step  5  .  Update  cost 


cost  =  cost  ♦  d  ♦  c 
k  i  i  ik 

If  ABRVT  <  1  ,  then  cost  =  cost  +  lp 
k  k  k  k  k 

If  ARRVT  >  u  ,  then  cost  =  cost  *  up 
k  k  k  k  k 


Step  6  .  let  i  =  k. 

If  N'  =  N,  then  go  to  next  step. 

Otherwise,  go  to  step  2. 

Step  7  .  Apply  the  Modified  Oropt  procedure  to  the 

current  tour.  Stop  when  no  further  improve¬ 
ments  can  be  found. 

End  of  algorithm  SCCO 

This  procedure  is  also  similar  to  the  cheapest 
selection,  cheapest  insertion  method  fcr  the  unconstrained 
TSP ,  except  that  the  nodes  with  time  windows  are  treated 
differently  from  the  time  free  nodes.  The  nodes  with  time 


windows  are  inserted  in  order  of  increasing  upper  time 
window  bounds.  The  time  free  nodes  are  inserted  between 
those  nodes  by  cheapest  selection  and  cheapest  insertion, 
for  as  long  as  the  upper  bound  of  the  time  windows  will 
allow. 

In  the  end,  a  Modified  Oropt  is  used  to  improve  the 
solution.  This  procedure  consider  only  those  exchanges  that 
would  result  in  a  node  being  inserted  netween  two  other 
nodes  in  the  current  tour. 

3.  SCAO 

This  algorithm  is  also  designed  for  the  time  window 
set  which  is  composed  of  some  tight  time  windows  and  some 
time  free  nodes. 

Algorithm  ;  SC  AO 

Input  :  Number  of  nodes,  x  and  y  co-ordinates  of  all 

nodes,  time  windows  for  all  nodes. 

output  :  Crdered  list  of  tour,  total  cost. 

Step  1  .  (initialization) 

Start  at  the  depot. 

Let  i=depot,  N *  =  (i),  cost  =  0. 

i 

Step  2.  Set  k  =  argmin  u  s.t.  j€N-N'. 

j 

If  k  is  time  free  node,  set  k  =  depot. 

Step  3  .  Insert  node  k  in  the  subtour  N'. 

Compute  AfiRVT 

k 

ARR  VT  =  ABfiVT  ♦  d  +  C 

k  i  i  ik 

Step  4  .  Insert  time  free  node  j  €  N-N*  between  nodes 

i  and  k  by  cheapest  insertion  and  greatest 

angle  (same  as  CCAO)  until  AfiRVT  does  not 

k 

exceed  u  . 

k 

t>1 


Step  5  .  Update  cost 


k 


cost  =  cost  «■  d  +  c 
k  i  i  ik 

If  AEBVT  <1  ,  then  cost  =  cost  ♦  lp 

k  k  k  k  k 

If  ARRVT  >  u  ,  then  cost  =  cost  ♦  up 
k  k  k  k  k 

Step  6  -  let  i  =  k. 

If  N'  =  N,  then  go  to  next  stop. 

Otherwise,  go  to  step  2. 

Step  7  .  Apply  the  Modified-Oropt  procedure  to  the 

current  tour.  Stop  when  no  further  improve¬ 
ments  can  be  fcund. 

End  of  algor  it  Jim  SCAO 

This  algorithm  is  same  as  SCCO  except  a  greatest 
angle  selection  in  stead  of  a  cheapest  selection  in  SCCC. 

C.  EXACT  SOIOTICV  TECHNIQUES  FOB  SOFT  TIME  UNDOES 
1.  State-Space  Relaxation  Procedure 

In  this  section  we  describe  a  state  space  relaxation 
procedure,  which  is  adapted  from  Christofides  et  al.  [Ref. 
5],  for  soft  time  windows.  They  only  considered  the  TSP  with 
hard  time  windows  and  without  time  windows.  The  dxfrerences 
are  as  follows.  Tne  waiting  cost  is  replaced  by  a  penalty 
cost  to  be  paid  in  the  early  arrival  case.  Late  arrival  is 
allowed,  but  a  penalty  cost  has  to  be  paid.  So  we  have  to 
calculate  tne  duration  and  the  penalty  cost  on  each  possible 
path  to  decide  the  least  cost  path  in  eacn  stage.  Re  denote 
the  penalty  cost  on  each  possible  path  as  PC  in  this 
section . 

Consider  the  TSP  defined  on  the  graph  G  =  {N ,  A)  with 
soft  time  window  constraints.  Let  S’  be  a  set  of  all  nodes 
except  starting  node.  Let  S  be  a  subset  of  S’.  Let  f  (S,  j) 


be  the  cost  of  the  least  cost  path  starting  at  node  1 
passing  through  every  node  of  S  and  finishing  at  node  j. 
Let  T(S,j)  he  the  total  duration  of  a  path  corresponding  to 
f(S,j).  Let  p(S,j)  be  the  predecessor  of  j  on  the  path 
corresponding  to  f  (S,  j) .  Let  lp(t)  be  the  early  arrival 
penalty  cost  function  and  up(t)  be  the  late  arrival  penalty 
cost  function.  For  a  given  S  and  j,  total  duration  of  a 
path  can  be  calculated  as 

T<S,j)  =  [  T(S-j,i)  ♦  d.  ♦  c.  .  ].  (4.1) 

i  il 

where  p  <S,  j)  =  i. 

In  equation  (4.1)  total  duration  of  the  least  cost  path 
passing  through  the  nodes  in  the  set  S  and  ending  in  node  j 
can  be  described  as  the  sum  of  three  terms:  the  first  is 
total  duration  of  the  least  cost  path  passing  through  the 
nodes  in  the  set  S-{j]  and  ending  in  node  i,  the  second  is 
the  time  required  to  spend  in  node  i,  and  the  third  is  the 
travel  time  frca  node  i  to  node  j.  The  dynamic  programming 
recursion  to  determine  the  least  cost  path  may  then  be 
stated  as 

f(S,j)  =  min  [  f(S-j,i)  «■  d  ♦  c  ♦  PC  ]  (4. 2) 

.  ies-j  i  ij 


where  T1 

=  [  l(S-j,i)  *  d. 

♦  c 

]• 

l 

ij 

PC 

=  0  ,  if  1  < 

T  1  < 

u 

lj 

=  lp  (1  -11)  ,  if 

■S 

T 1  J 

=  up  (TT-u  )  ,  if 

initialization: 

I  1  > 

3 

u 

3 

with  the 

£  <{j}  r  j)  =  c  ,  if 

1  < 

c 

< 

u 

13  ,  „ 

j 

1  j 

j 

=  c  ♦  lp  (1 

"c,  ) 

9 

if 

c 

1  j 

3 

if 

lj 

=  c  ♦  up  (c 

.  .~u  > 

9 

c 

1  j 

1  3  3 

- 

1  j 

Finally,  the  optimum  solution  can  be  calculated  as 


ain  [  f(S*,i)  ♦  d 


♦  c .  ]. 


Since  the  computer  storage  requirements  are 
increased  exponentially  with  the  size  of  the  problem,  this 
method  is  limited  to  small  problems.  For  relaxing  this 
limitation,  a  state  space  relaxation  procedure  can  be  used 
same  as  Chapter  III. 

Consider  the  dynamic  programming  formulation  (4.2) 
The  state  variable  in  that  formulation  is  (S,j),  and  the 
stage  is  the  cardinality  of  S.  Let  g  (S)  be  a  mapping  from 
the  domain  of  (S,j)  to  some  other  vector  space  (g(S),g). 
Let : 


H(g(S),j)  =  {  (g(S-j),i)  |  i  €  s-  j  } 


(4.3) 


Since  we  are  interested  in  lower  bounds  to  the  TSP 
with  time  constraints,  H(g(S),j)  in  (4.3)  may  be  replaced  by 
any  larger  set  that  is  easier  to  compute.  Thus,  H(g(S),j) 
can  be  defined  by  the  following  equation: 

H(g(S),j)  =  {  (g(S-j),i)  I  i£E(g(S),j)  }  (4.4) 

where  S-  j  £  E  (g  (S)  ,  j)  . 

For  calculating  the  lower  bound  of  the  problem,  recursion 
(4.  1)  can  be  changed  to  the  following  equation: 


.  T  (g(S)  ,  j)  =  [  T  (g  (S- j)  ,i)  ♦  d .  ♦  c.  .  ] 

i  13 

where  p(g(S)  , j)  =  i. 

Recursion  (4.2)  may  be  stated  as 


(4.5) 


f(<J(S),j)  *  min  [f  (g  (S- j)  ,i)  ♦  d  +  c.  ♦  PC]  (4.6) 

(9(S-3) ,i) e  H  (g  (  5)  ,  g)  r  ig 

=  min  (f  (g  (S- j)  , i)  ♦  d  ♦  c  *  PC]  (4.7) 

i  £  E(g(S),j)  i  ig 

where  T 1  =  £  7(g(S-j)  ,i)  ♦  d  ♦  c  ], 

i  ij 

PC=0  ,  if  1  <  T  1  <  u  . 

=  lp  (1  -T 1 )  ,  if  T  1  <  1D 


with  the  initialization: 


f  (9(j)  #  j)  =  c 

=  c  *■ 


,  if  1  <  c  <  u 

,  .  lpU  -c  .?  ,  if3  c  j<  1  . 

id  A  .  ]  i]  ,  ..  1 i  .  i 

=C  ♦  Up(C  -U  )  ,  if  c  >u 

Ij  Ij  j  1j  j 


Finally,  the  optimum  solution  can  be  calculated  as 

r  f  (g  (s • )  #i )  ♦  d. 

» 1)  J 


am 

i  e  E  (g  (  S) 


♦  c  ]. 
i  1 


The  mapping  can  be  selected  frca  any  separable  func¬ 
tion.  He  used  a  mapping  function  (3.7),  which  is  proposed  by 
Christofides  et  al. ,  same  as  Chapter  III.  Then  equation 
(4.  5)  becomes  : 


T  (|S|  ,  j)  =  [  T  ( |  S|  -  1,i)  ♦  d.  ♦  c  ]. 

i  ij 


(4.8) 


where  p(|S  |  ,  j)  =  i 

fiecursicn  (4.7)  may  be  stated  as: 


f  (|S|  ,  j)  =.  min  [  f  (<S|-1 
i  e  E(1S(  ,  g) 


-1,i)*  d  ♦  c  +  PC  ] 
i  ij 


(4.9) 


where  T1  =  [  T(|S|-1,i)  +  d  ♦  c  ], 

i  ij 


PC  =  0 


0  ,  if 

1  _  < 

T  1 

lp  (1  -T1) 

»  if 

T  1 

up  (T?-U  ) 

#  if 

T  1 

>  u 


with  the  initialization: 
f  <1,  j)  =  c 


=  c 


Ij 


<  u 


,  if  1  <  c 

j  1  j  j 

♦  lp(l  -C  )  ,  if  C  .  <  1. 

1]  J  1]  1 g  j 

=  c  +  up (c  -u  )  ,  if  c  >  u 
Ij  Ij  j  Ij  j 


Finally,  the  optimum  solution  can  be  calculated  as 


min  [  f  { |  S'  |  ,i)  ♦  d  +  c  ], 

i  e  E ( |  N |,1)  i  il 


2 .  Additional  Condi tion 


In  the  previous  section,  we  discussed  a  state  space 
relaxation  procedure  which  is  adapted  from  christofides  et 
al.fBef.  5].  That  procedure  provides  a  lower  bound  on  the 
ISP  with  soft  tine  window  constraints.  The  additional 
condition  to  avoid  loops  formed  by  three  consecutive  nodes 
was  used  to  get  a  better  bound  [Eef.  5],  This  can  be  done 
in  the  following  way. 

Let  k  =  |  S  | .  Let  f(k,j,1)  be  the  cost  of  tne  least 
cost  path  frcm  the  initial  state  to  state  (k,  j)  without 
loops  formed  by  three  consecutive  nodes.  Let  f(k,j,2)  be 
the  cost  of  the  second  least  cost  path  from  the  initial 
state  tc  state  (k,  j)  without  loops  formed  by  three  consecu¬ 
tive  nodes.  Let  p(k,j,m)  be  the  predecessor  of  j  on  the 
path  corresponding  to  f{k,j,m).  With  the  above  definition, 
equation  (4.8)  becomes: 

T(krj,^l,)  =  [  T(k-1,i,m)  +  d  +  c  ],  m'  =  1,2  (4.10) 

i  ij 


where  F  (k,  j  ,  n* )  =  i 


m  =  1  ,  if  p(k-  1,i,1)  *  j 
=  2  ,  otherwise. 


With  the  initialization: 
T  (1.  j,  D  =  c  . 


and 


T(1,j,2)  =  »  . 

Recursion  for  f(k,j,1)  can  be  calculated  in  the  following 
way.  Let: 

T’  (k,j,m)  =  £T(K-1,i,m)  +  d  ♦  c  ],  m  =  1,2 

i  ij 


This  gives  us: 


f(k,  ],1)  =  min  ff(k-1,i,m)  ♦  d  +  c  ♦  PC  ]  (4.11) 

i  e  E(k,  j)  i  ij 


where  PC  =  0  ,  if  1  <  T*  (k,j,m)  <  a 

=  lp  (1  -T  •  (k  ,  j3m ) ) ,  if  T '  (k,  j,  mj'  <  1 


=  up  <T  •*  (k,  j  ,  m)~-u  )  ,  if  T'(k,j,m)  >  u3 


m  =  1  ,  if  p  (k-1,i,1) 
=  2  ,  otherwise. 


With  the  initialization: 
f  (1,i,  1)  =  c 


1i 


,  if  1  <  c  <  u 

i  1  i  i 

=  c  ♦  lp(l .  -c  )  ,  if  c 

1i  i  1i  1i 

=  c  ♦  up(c  -a.)  ,  if  c 

1i  1i  i  1i 


(4-12) 


<  1 


>  u 


Recursion  for  f(k,  j,2)  can  be  written  in  the  following  way: 


i  e  E  (k, 

i*P(k, j,  1) 


where  PC  =  0 


«■  d  ♦ 

c 

i 

i  j 

)  *  u 

(k, 

<  1 

(k,  j,  m) 

>  u 

a  =  1  „  if  p  (k-  1,i  ,  1)  *  j 
=  2  ,  otherwise. 


With  the  initialization: 
f  ( 1,  i,  2)  =  « 

Finally,  the  optimum  solution  can  be  calculated  as 


(4.14) 


ain  [  f  ( |S»  | ,  i,  1 )  ♦  d  +  c  ]. 
i  e  E  ( J  N  |  ,  1)  i  i  1  J 


(4.15) 


Since  the  additional  condition  can  avoid  considera¬ 
tion  of  a  useful  lower  bound,  we  considered  f(k-1,i,2)  in 
recursion  (4.11)  and  (4.13)  only  when  the  predecessor  of  i 
on  the  path  ccrres pondinq  to  f(k-1,i,1)  is  j.  If  we  do  not 
consider  the  second  least  cost  path  in  case  of  p  ( k- 1  ,i  ,  1)  =  j  , 
then  f(g(S)  ,j)  does  not  guarantee  the  lower  bound  of  f  (S,j) 


For  this  example,  let's  consider  4  node  TSP  with 
time  constraints.  Node  A  is  the  starting  node.  D  is  the 
time  free  node.  The  lower  bound  of  node  B  is  9,  the  upper 
bound  of  node  B  is  11,  the  lower  bound  of  node  C  is  19,  and 
the  upper  bound  of  node  C  is  21.  Suppose  service  time  at 
each  node  is  2ero,  lp(t)=t,  and  up(t)  =  5t.  Figure  3.7  shows 
an  optimal  route  for  this  problem.  Frcm  equation  (4.10)  and 
(4.  12)  we  can  get: 

f(1,B,1)  =  10,  T(1,B,1)  =  10,  p  (  1,  B,  1 )  =  A; 

f  (1#C,  1  )  =  1  9#  T  ( 1 ,  C,  1  )  =  14.14,  p  ( 1  ,C,  1)  =  A; 

f  ( 1,  D,  1 )  =  7.07,  T  (  1,D  ,  1)  =  7.07,  p(1,D,1)  =  A. 

New  applying  equation  (4.10)  and  (4.11)  recursively  with 
i= 1 ,  for  k=  2  we  can  get: 

f  (2,3,1)  =  min  [  94.7,  29.84  ]  =  29.84, 
i  e  {C  ,  D) 

T  (2,  B  ,  1 )  =  14.14,  p  (2,  B,  1)  =  D. 


Simila  rly , 


f(2,C,1)  =  19,  T  (2,  C,  1  )  =  14.14,  p  (2,C,  1)  =  D; 
f  (2, 0,1)  =  17.07,  1  (2,  D,  1)  =  17.07,  p  (2,0,1)  -  B. 


Fcr  k  =  3, 


f  (3,  B,1)  =.mincj  94.7  ]  =  94.  7, 


T  (3,  B  ,  1 )  =  24.14,  p  (3 ,  B  ,  1)  =  C. 

Similarly, 

f  (3,C ,  1 )  =  39.84,  T(3,C,1)  =  24.14,  p(3,C,1)  =  D; 
f  (3,D,  1)  =  «>  . 

We  can  see  easily  that  f(3,D,1)  is  not  a  lower  bound  of 
f  <{B,C,C)  ,D  )  . 


3.  Branch  and  Boand  Procedure 

He  used  the  same  branch  and  bound  procedure  used  to 
eliminate  subtours  in  the  solution  of  the  state  space  relax¬ 
ation  procedure  in  Chapter  III.C.  3.  The  SCCO  heuristic, 
which  was  described  in  section  B.2,  was  used  as  an  initial 
upper  bound,  and  the  lower  bound  was  obtained  from  equation 
(4-15)  . 

We  present  the  results  of  our  computational  experi¬ 
ence  with  the  algorithms  of  this  Chapter  in  the  next 
Chapter . 


V.  COHPOT ATI  ORAL  El PEBIENCB 


A.  TEST  PROBLEMS 

Four  sets  of  test  data  are  used  in  this  thesis.  Test 
problem  number  [ 1 ]  is  taken  from  Sedgewick  [fief. 19:  p.309]. 
The  other  problems,  numbered  [2],  [3]  and  £4],  are  from 
Appendix  9.1  cf  Eilon  et  al.  *  s  text  ffief.  21].  These  test 
problems  are  shown  in  Appendices  A,B,C  respectively.  These 
published  problems  contain  node  and  depot  locations,  but 
they  do  not  include  time  windows. 

We  constructed  time  windows  for  test  problems 
[1],£2],£3]  by  first  using  the  CCAO  heuristic  on  the  uncon¬ 
strained  ISP.  Time  windows  were  then  placed  acout  each  node 
such  that  the  CCAO  route  was  feasible.  The  idea  for  gener¬ 
ating  time  windows  in  this  way  comes  from  Baker  ffief.  25], 
who  used  the  unconstrained  Nearest  Neighbor  heuristic  as  his 
starting  point  instead  of  CCAO. 

The  time  window  widths  were  set  to  varying  sizes  ranging 
from  3  to  14.  Seme  of  the  time  windows  were  rairly  tight 
while  others  overlapped.  This  is  in  contrast  to  Baker's 
work,  where  all  the  time  windows  have  width  equal  2  units. 

The  last  problem  number  £4]  is  the  same  as  test  problem 
£3],  except  that  the  time  windows  were  constructed  from  a 
Nearest  Neighbor  solution  to  the  unconstrained  traveling 
salesman  problem,  as  in  Baker  [Ref.  25].  Figure  5.1 
displays  the  CCAO  solution  for  test  problem  [3]  and  Figure 
5.2  illustrates  the  unconstrained  Nearest  Neighbor  solution 
for  the  test  problem  £4].  Me  found  a  small  error  in  Baker's 
TSP  solution  for  the  Nearest  Neighbor  [Ret.  25],  in  that  the 
nearest  node  from  node  16  is  node  17,  not  node  13.  The 
resulting  cost  is  actually  higher,  it  is  312.09,  not  310.22. 


Figure  5.1  Unconstrained  Solution  Obtained  by  CCAO. 

Each  of  the  four  sets  of  test  data  was  used  to  create 
four  test  problems.  The  separate  instances  differed  in  the 
percentages  of  time  window  constraints  that  were  chosen  to 
be  in  effect-  The  four  cases  were  100%,  90%,  75%,  and  50%. 
We  refer  to  this  percentage  as  the  "time-window  percentage”. 

A  random  number  generator  was  used  to  decide  which  nodes 
would  have  time  windows.  Test  problems  for  the  time  window 
constrained  TSP  are  shown  in  Appendices  G  through  V.  The 


Figure  5.2  ancon strained  Solution  Obtained  by 
Nearest  Neigh tor  Heuristic. 

Fenalty  cost  factors  can  be  varied  depending  on  real  world 
problems.  We  used  2  and  5  as  the  lower  and  upper  penalty 
cost  factor.  Also  we  set  the  service  time  at  each  node  to  0 
to  make  it  easy  to  construct  the  time  windows. 

The  computational  results  are  presented  in  T.ables  II  and 
III.  The  figures  reported  represent  results  of  our  test 
runs  for  each  test  case. 


B.  COfl POTATION AX  BESOXIS 


1.  Hard  lime  Windows 

As  noted  in  section  3,  the  Nearest  Neighbor 
heuristic  often  cannot  solve  the  problem,  because  the 
arrival  time  of  the  nearest  node  frequently  violates  the 
upper  bound.  However  in  test  problem  [4],  the  results  of 
the  Nearest  Neighbor  are  the  same  as  in  the  unconstrained 
problem,  because  this  problem  itself  was  constructed  by  a 
Nearest  Neighbor  heuristic. 

Generally,  the  SCCO  and  SCAO  heuristic  can  be  easily 
applied  to  the  hard  time  window  TSP.  According  to  our 
experiments,  if  the  time  window  width  becomes  large  relative 
to  the  travel  time  between  nodes  on  the  optimal  uncon¬ 
strained  TSP  route,  then  the  lower  percentage  time  window 
problems  become  more  difficult  to  satisfy.  This  phenomenon 
can  be  seen  in  test  problem  [1].  That  is  because  the  other 
nodes  in  the  cptimal  route  for  the  unconstrained  TSP  problem 
could  be  inserted  without  causing  violation  of  the  upper 
bound. 

The  STACK  heuristic  takes  slightly  more  time  than 
the  other  heuristics.  It  achieved  lcwer  accuracy  in  test 
problems  [1  ]  and  [3]  in  the  50  percent  time  window  case. 
The  exact  algorithm  can  find  the  exact  answer  in  most  prob¬ 
lems,  but  when  there  are  fewer  windows  in  effect,  it  takes 
mere  computation  time.  It  cannot  solve  the  50  percent  time 
window  problems  [2]  and  [4]  within  180  seconds. 


Number  lime  Nearest  SCCG  SCAD  SLACK  Exact 

Problem  of  Windows  Neighbor 
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2.  Soft  lime  Windows 

All  of  the  methods  tested  for  soft  time  windows  were 
able  to  find  some  answer  to  every  problem  within  reasonable 
computing  time,  except  for  two  instances  with  the  exact 
algorithm.  With  the  Nearest  Neighbor  heuristic,  the  quality 
of  solution  is  not  desirable.  In  general,  the  lower  time 
window  percentage  problems  have  lower  solution  quality.  As 
in  the  hard  window  case,  on  test  problem  [4],  the  results  of 
the  Nearest  Neighbor  heuristic  coincide  with  the  uncon- 
strained  TSP  heuristic,  because  this  problem  itself  was 
constructed  b;  a  Nearest  Neighbor. 

As  in  the  hard  time  window  problem,  SCCO  and  SCAC 
generally  find  an  optimal  solution  except  for  one  problem 
witn  50  percent  time  windows.  In  test  problem  £1]  with  50 
percent  time  windows,  the  SCCO  and  SCAO  values  were  215.686, 
165.544  respectively.  The  exact  algorithm  could  not  solve 
the  two  test  problems  with  50  percent  time  windows  within 
180  seconds.  The  reason  is  that  the  solutions  of  the  state 
space  relaxations  have  many  subtours  and  it  takes  a  long 
time  to  eliminate  these  subtours. 

With  both  hard  and  soft  time  windows,  the  results 
are  sensitive  tc  the  percentage,  width  and  position  of  the 
time  windows.  In  most  problems,  the  fewer  time  windows 
there  are,  the  lower  the  accuracy  of  the  heuristics. 


.  Number  lime  Nearest  SCCO  SCAO  Exact 

Problem  of  windows  Neighbor 
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VI.  COICLOSXOMS  AMD  B SCO MM EM PATIO MS 


This  thesis  has  preseated  some  heuristics  and  exact 
algorithms  for  the  solution  cf  traveling  salesman  problem 
with  ti se  window  constraints.  Me  considered  two  different 
kinds  of  time  window  constraints  :  hard  time  windows  and 

soft  time  windows.  Hard  time  windows  are  inviolable, 
whereas  soft  windows  may  be  violated  at  a  cost. 

For  both  hard  time  windows  and  soft  time  windows,  we 
developed  some  new  heuristics,  SCCC  and  SCAO,  which  are 
modifications  of  Stewart's  unconstrained  TSP  heuristics 
[Ref.  6]  CCCO  and  CCAO.  Also  for  the  hard  time  window  only, 
we  developed  the  SLACK  heuristic.  Me  also  developed  an 
exact  algorithm  for  both  hard  and  soft  window  using  state 
space  relaxation  dynamic  programming  and  branch  and  bound  as 
proposed  by  Christofides  et  al.  [Ref.  5]. 

The  procedures  were  shown  to  be  effective  on  some  moder¬ 
ately  small  sized  problems.  A  Nearest  Neighbor  heuristic 
was  also  developed,  but  it  was  often  unable  to  solve  the 
problem  with  hard  time  windows,  and  it  found  very  low 
quality  solutions  with  soft  time  windows.  This  experience 
is  consistent  with  the  findings  of  ethers  [Ref.  7]  who 
determined  that  the  Nearest  Neighbor  heuristic  does  not 
perform  well  cn  the  unconstrained  TSP. 

The  SCCQ  and  SCAO  are  generally  effective  cn  most  of  the 
small  sized  problems  we  tested,  except  for  the  problems  in 
which  less  than  half  the  nodes  have  time  windows.  Furtner 
research  is  needed  in  order  to  satisfactorily  solve  these 
problems.  Another  problem  difficulty  that  may  require  more 
research  is  dealing  with  wider  time  windows. 

The  SLACK  heuristic  which  is  used  only  with  hard  time 
windows  is  slightly  slower  than  the  other  heuristics. 
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Particularly,  in  the  lower  time  window  percentage  problems, 
the  accuracy  becomes  lower. 

The  exact  algorithm  succeeded  in  solving  14  of  the  16 
test  problems  to  optimality,  but  it  was  too  slow  to  use  in 
most  of  the  lower  time  window  percentage  problems.  This 
algorithm's  performance  also  depends  upon  the  quality  of  the 
upper  bound  which  is  obtained  from  the  heuristic. 
Additional  research  is  needed  to  reduce  computation  time, 
but  a  working  program  for  at  least  some  problems  has 
resulted  from  this  effort. 
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Depot  co-ordinates  :  (326,1  81) 

problem  source  :  Ellon  et  al.  [6ef.  21] 
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APPEBDII  B 
TEST  PBOBLEB  [6] 
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Depot  co-ordinates  :  (40,40) 

problea  source  :  Ellon  et  al.  [ fief  .21] 


APPBMDII  F 

TEST  PRO BLEU  FOB  THE  SCCO 
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Depot  co-ordinates  :  <12,10) 

CL  =  2. 0,  CO  =5.0 

problem  source 

node  locations  :  Sedgewick  [Ref.  19] 
time  windows  :  see  Chapter  V. 


|  node 

X 

y 

time 

l(i) 

window 

u(i) 

node 

X 

y 

time 

l(i) 

window 

a  (i) 

t  1 

3 

9 

25 

32 

11 

10 

13 

10 

17 

1  -2 

11 

1 

46 

53 

12 

16 

14 

2 

9 

1  3 

6 

8 

27 

36 

13 

15 

2 

51 

59 

1  4 

4 

3 

- 

- 

14 

13 

16 

5 

13 

1  5 

5 

15 

18 

28 

15 

2 

12 

22 

30 

1  6 

8 

11 

14 

23 

16 

12 

10 

- 

- 

|  7 

6 

4 

35 

43 

1  3 

7 

4 

42 

49 

1  9 

9 

7 

58 

68 

1  io 

14 

5 

53 

64 

Depot  co 

-ordinates  : 

(12,10) 

CL  =  2.0 

,  C(J  =5.0 

problem 

source 

node 

locations  : 

Sedgewick  [Ref 

time 

windows  : 

see  Chapter  V. 

APPE1DII  I 
TEST  PBOBIBH  [1-3] 
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CL  =  2.  0,  CU  =5.0 
problem  source 

node  locations  :  Sedgewick  [Ref 
tine  windows  :  see  Chapter  V. 
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APP  B1DII  J 
TEST  PfiOELBH  £1-4] 
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:  see  Chapter  V. 
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APPE1DII  a 
TEST  PROBLEM  [2-3] 
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Depot  co-ordinates  ;  (326,181) 

CL  =  2.  0,  CU  =  5.0 

problem  source 

node  locations  :  Eilon  et  al.  [fief.  21] 
time  windows  :  see  Chapter  V. 
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Depot  co-ordinates  :  (143,215) 

CL  =  2.0  ,  CD  =  5-0 

problem  source 

node  locations  :  Eilon  et  al  [Kef.  21]. 
tine  windows  :  see  Chapter  V. 
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APPEHDII  2 
TEST  PROBLEM  [3-3] 
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APPEIDII  a 
TEST  PRO  BLEU  £3-4] 
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Depot  co-ordinates  :  (145,215) 

CL  =  2.0  ,  CO  =  5.0 

proJblcm  source 

node  locations  :  Eilon  et  al.  [Ref.  21]. 
time  windows  :  see  Chapter  V. 
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1PPE1DII  S 
ZEST  PROBLEM  £4-1] 
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Depot  co-ordinates  :  (145,215) 

CL  =  2.0,  CO  =  5.0 
problem  soarce 

node  locations  :  Eilon  et  al  [fief.  21] 
tine  windows  :  see  Chapter  V. 
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tine  windows  :  see  Chapter  V. 
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